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Abstract: This paper, which expands on a talk given at the International Workshop on Fat 
Points, February 9-12, 2000, in Naples, Italy, surveys problems and progress on certain problems 
involving numerical characters for ideals I{Z) defining fat points subschemes Z = m\pi + ■ ■ ■ + 
m n p n C P 2 for general points pi. In addition to presenting some new results, a collection 
of MACAULAY 2 scripts for computing actual or conjectured values of (or bounds on) these 
characters is included. One such script, f indres, for example, computes the syzygy modules in 
a minimal free resolution of the ideal I(Z) for any such Z with n < 8; since f indres does not 
rely on a Grobner basis calculation, it is much faster than routines that do. 



I. Introduction 

This paper surveys work on certain problems involving fat points subschemes of P 2 . 
To encourage experimentation, I have included a number of MACAULAY 2 scripts for 
doing explicit calculations. To simplify using them, I've included them in the TjrjXfile for 
this paper in a verbatim listing, without any intervening IgX control sequences. Thus if 
you have (or obtain, from, say, [http : / / www . math . unl . edu/ ^bharbour /Survey . tex|) the 



Tf^Xlisting for this paper, you can simply copy the lines for the necessary MACAULAY 2 
scripts from this paper directly into MACAULAY, without any additional editing. 

Although the most general definition of a fat points subscheme involves the notion of 
infinitely near points (see [H6]), it is simpler here to define a fat points subscheme of P 2 
to be a subscheme Z defined by a homogeneous ideal / C R of the form I(pi) mi PI • • • H 
J(p n ) m ™, where pi, . . . ,p n are distinct points of P 2 , mi, . . . , m n are nonnegative integers 
and R = k[P 2 ] is the homogeneous coordinate ring of P 2 (i.e., a polynomial ring in 3 
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variables, x, y and z, over an algebraically closed field k). It is convenient to denote Z by 
Z = m\p\ + • ■ ■ + m n p n and to denote / by I(Z). 

For another perspective, a homogeneous polynomial / G R is in I(Z) if and only if 
mult Pi (/) > mi for all i, where mult Pi (/) denotes the multiplicity of / at pi (this being the 
least t such that if l\ and I2 are linear forms defining lines which meet at pi and nowhere 
else then / is in the tth power (Zi, I2Y of the ideal (Zi, Z2) G R). 

It is important for what follows to note that I(Z) is a homogeneous ideal, hence I(Z) 
is the direct sum of its homogeneous components I(Z) t = R t C\I(Z) (where for each integer 
t, R t denotes the k- vector space span of all homogeneous polynomials of R of degree t). 

1.1. Numerical Characters 

The work I am interested in here concerns certain numerical characters of ideals 
I(mipi + • ■ • + m n p n ) C R which take a constant value on some nonempty open sub- 
set of points (pi, . . . ,p n ) G (P 2 ) n . Thus we will usually consider fat points subschemes 
Z = m\p\ + ■ ■ ■ + m n p n for which the points pi G P 2 are general. (Saying that something 
is true for Z — ui\p\ + ■ ■ ■ + m n p n for general points pi, is the same as saying that it holds 
for some open subset of points (pi, ... , p n ) G (P 2 ) n .) In order to establish a result for gen- 
eral points, one typically establishes it for some particular special choice of the points and 
then argues by semicontinuity. (To justify using semicontinuity, even for specializations 
to infinitely near points, see my 1982 thesis, the relevant parts of which were published 
in [H2]; alternatively, for specializations keeping the points distinct, see [P].) Thus we will 
sometimes consider situations for which the points pi are in some special position. 

Given a fat points subscheme Z = m\p\ + • • ■ + m n p n and its ideal / = I(Z), among 
the numerical characters which have seen attention by various researchers over the years 
are the following: 

• ot(Z), the least degree t such that I(Z) t 7^ 0; 

• P(Z), the least degree t such that the zero locus of I(Z)t is zero dimensional; 

• hz, the Hilbert function of I(Z) (i.e., the function whose value hz(t) for each 
degree t is the k- vector space dimension of I(Z) t ); 

• t(Z), the least degree t > such that hz(t) = Pz(t), where Pz is the Hilbert 
polynomial of Z (which is simply Pz(s) = (s 2 + 3s + 2 — £\ mj(mj + l))/2); 

• v t (Z), the number of generators of I(Z) in degree t in any minimal set of homo- 
geneous generators. 

(In cases where it is understood which Z is meant, I will sometimes write a or /?, etc., for 
the more explicit but more cumbersome a(Z), etc.) 

The most fundamental characters are a, h and v t . For example, hz immediately 
determines oe(Z) and t(Z). Moreover, if one can compute hz for any Z then one can also 
determine (3 for any particular Z. (This is because t < (3(Z) if and only if either t < ot(Z), 
or t > a(Z) and there exists some nonzero Y = m' x p\-\-- • ■ + m' n p n with < vn! i < mi for all 
i such that hz(t) = hz-yit — a(Y)). The idea is that for a(Z) < t < (3(Z), every element 
of I(Z)t is divisible by nonconstant homogeneous polynomials / which define divisors in 
the fixed locus of the linear system I(Z) t . Any such / spans I(Y) a (Y) for an appropriate 
Y, as above. Since there are only finitely many vn! x p\ + ■ ■ ■ + m' n p n with < m\ < mi, 
one can in principle check whether any such Y exists, as long as one can always compute 
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h.) Similarly, if one can compute a(Z) for any Z then one can also determine hz for any 
particular Z . (Here's how: to compute hz(t) for some t and some Z = m\p\ + ■ ■ ■ + m n p n , 
let Z = Z and for each j > let Zj = Z + qi + ■ ■ ■ + qj, where q±, . . . , qj are general. 
Since each additional point qi imposes one additional condition on forms of degree t up to 
the point where no forms remain, we see that hz(t) = i where i is the least j such that 
a(Zj) > t.) 

Knowing a(Z) for a particular Z sometimes also means we know vt{Z") for all t. 
Indeed, conjectures about the values of a, and, for certain Z, of u t , are made below (see 
Conjecture II.3.1). There are examples of Z for which, if a(Z) is what it is conjectured to 
be, then so are all v t {Z}. (For example, take Z = m(pi + ■ ■ ■ + p n ) with general points Pi, 
where n > 9 is an even square and m is sufficiently large; see Example 5.2 and Theorem 
2.5, both of [HHF].) Nonetheless, knowing a in general does not seem to be enough to 
determine v t for all t, but if one knows v t for all t, then one can always compute a and 
hence all of the other characters listed. (This is because the least t such that vt(Z) > 
is t = a(Z); moreover, z/ a (^)(Z) = hz(a(Z)).) Thus the characters u t are perhaps even 
more fundamental than the other characters discussed above. 

The characters vt{Z") are also interesting due to their connection to minimal free 
graded resolutions of I(Z). A minimal free graded resolution of I(Z) is an exact sequence 
— > Fi(Z) — > Fq — > I(Z) — > in which Fi(Z) and Fq(Z) are free graded R- modules. It 
turns out, up to isomorphism as graded -R-modules, that Fo is (BtR[— t] Ut ^ and F\(Z) 
is ® t R[-t] St ^ Z \ where the characters s t {Z) are defined via v t {Z) - s t {Z) = A 3 hz(t) 
[FHH]. Here A denotes the difference operator (so for any function / : Z — > Z, we have 
= f(t) — f(t — 1)), and R[iy denotes the direct sum of j copies of the module R 
itself but taken with the grading defined by R[i]t = Rt+i- 

1.2. Connection to Geometry 

Additional interest in these characters (and essential techniques in studying them) 
comes from their connections to geometry. Given distinct points pi,...,p n £ P 2 , let 
7T : X — > P 2 be the birational morphism obtained by blowing the points up. Thus n is the 
unique morphism where X is a smooth and irreducible rational surface such that, away 
from the points pi, it is an isomorphism and such that for each i, ir~ 1 {pi) is a smooth 
rational curve Ei. It is known that the divisor class group C1(X) is a free abelian group 
on the classes [Ei] of the divisors Ei and on the class [Eq], where Eq = 7r _1 (L), L being 
any line in P 2 not passing through any of the points Pi. Thus for any divisor D on X we 
have [D] = Yl7=o a i[^i] f° r some integers a^. 

It will be useful later to recall the intersection form on C\(X). This is a symmetric 
bilinear form denoted for elements [C] and [D] of C\(X) by [C] ■ [D] and determined by 
requiring that [Ei] ■ [Ej] is if i ^ j, 1 if i = j = and — 1 if i = j > 0. If C and D are 
curves on X such that C fl D is finite and transverse, then [C] ■ [D] = \C(~)D\; i.e., [C] ■ [D] 
is just the number of points of intersection of C with D. If C = D, it is convenient to 
denote [C] ■ [D] by [C] 2 . 

Now, for a divisor D, let Ox(D) denote the associated line bundle. Given a fat points 
subscheme Z = m\p\ + • • ■ + m n p n , it turns out for all t that hz(t) — h°(X, Ox(Ft(Z))), 
where F t (Z) is the divisor tEo — (m\Ei + ■ — \-m n E n ) and h°(X, Ox(F t (Z))) denotes the 
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dimension of the Oth cohomology group H°(X, Ox{F t (Z))) of the sheaf Ox(F t (Z)) (i.e., 
h°(X,Ox(Ft(Z))) is the dimension of the space of global sections of Ox(F t (Z))). Thus 
a(Z) is the least t such that h°(X, O x{F t (Z))) > and t(Z) is the least t > such that 
h°{X,O x {F t {Z))) = P z {t). Moreover, note that P z (t) = (F t (Z) 2 - K x • F t (Z))/2 + 1. 
By Riemann-Roch we have 

/i°(X, 0x(F t (Z))) - ^(X, Ox(Ft(^))) + h 2 (X, O x (F t (Z))) = P z (t), 

and by duality we know h 2 (X, O x {F t (Z))) = for t > -3, so h°(X, O x (F t (Z))) - 
h l {X,O x {F t {Z))) = P z (t) for all t > 0. Thus h z (t) = P z {t) + h\X, O x {F t {Z))) = 
dimRt - (Y /i m l (m t + l)/2 - h}{X,O x {F t {Z)))). Now, I{Z) t is precisely what is left 
from R t after imposing for each i the condition of vanishing at pi to order at least 
mf, what the previous equation is saying is that the number of conditions imposed is 
Ei ™>i( m i + i)/ 2 - h\X, O x (F t (Z))). For all t sufficiently large, h\X, O x (F t (Z))) = 
so a total of J2i m i( m i + conditions are imposed. For smaller t, /i^J, Cj(F t (Z))) 
measures the extent to which these Ei m i( m i + conditions fail to be independent, 
and we can regard r(Z) as the least degree in which the conditions imposed become inde- 
pendent. 

Likewise, the characters v t can be understood from two perspectives. There is a 
natural map m(Z) : I(Z) t ®fc Ri — > I(Z) t +i given by multiplication, and u t +i(Z) is 
just the dimension of the cokernel of the map ^t(Z). Corresponding to this map ^t{Z) 
we have in a natural way a map /j(F t (Z)) : H°(X,O x (F t (Z))) (g) k H (X,Ox(E )) -> 
H°(X, O x (F t+1 (Z))) : and indeed u t+1 (Z) = dimcok(^(F t (Z))). 

II. Resolutions 

Given a fat points subscheme Z = rriipi + ■ ■ ■ + m n p n , much current work concerns 
either computing or bounding one or another of the numerical characters cited above. 
Some of the oldest such work concerned bounding the characters v t . 

II. 1. Dubreil and Campanella Bounds 

Dubreil [Dub] obtained two bounds on the minimum number E» v i(Z) of homogeneous 
generators of I(Z): 

Theorem II. 1.1: Let Z = m\p\ + • • • + m n p n be a fat points subscheme of P 2 with 
distinct points p l . Then ^ v%{Z) < a(Z) + (3{Z) - r(Z) < a(Z) + 1. 

Sketch of proof: The inequality < a (Z) + 1 follows immediately from the 

Hilbert-Burch Theorem. Here is a more elementary proof. Given R = k[x,y,z], we may 
assume that x, y and z define general lines in P 2 (which, in particular, do not contain any 
of the points pi). It is then easy to see that the image J = xI(Z) t + yI(Z) t of the map 
xI{Z)t®yI{Z) t — > I(Z) t +i has dimension 2h z (t) —h z (t — l). Since J has a base point (all 
elements of J vanish at the common point of vanishing of x and y), we see for all t > oc(Z) 
that xI(Z) t + yI(Z) t cannot contain zI(Z) t . Hence for t > a(Z) the image of fJ^t(Z) has 
dimension at least 2h z (t) — h z (t — 1) + 1, hence ut+i(Z) = dimcok(/it(Z)) < h z (t + 1) — 
(2h z (t) - h z (t -!) + !) = A 2 h z (t + !)-!, while of course u t (Z) = h z (t) = A 2 h z (t) 



Survey of Fat Points on P 



5 



for t = a(Z). Summing for % = a(Z) to any N sufficiently large so that Vj(Z) = and 
h z (J) = p z{j) for j >N-1, we obtain < 1 + T,i(& 2 h z (t) - 1) = l + Ah z (N)- 

(N - (a(Z) - 1 )) = P Z (N) - P Z (N - 1) - N + a(Z) = N + 1 - N + a(Z) = a(Z) + 1. 

The foregoing proof is based on an argument given by Campanella [Cam]. Using a 
result of [GGR], Campanella there also gives a similar but slightly more refined bound, 
vt+i(Z) < A 2 h z (t + 1) - e t , where e t is for t < a(Z), 1 for a(Z) < t < (3{Z) and 2 
for P(Z) < t < t(Z). Summing these refined bounds for i from a(Z) to t(Z) + 1 gives 
i>i(Z) < a(Z) + (3(Z) — t(Z). (This argument requires that one knows that Vi(Z) = 
for i > t(Z) + 1, but this is true and well known; see [DGM]. It is also not hard to see 
this directly, at least from the point of view of the surface X obtained by blowing up the 
points pi. Let [E ], . . . , [E n ] be the corresponding basis for C1(X), as discussed in Section 
I. The statement we need to prove is then that n(Ft(Z)) is surjective if t > r(Z). But 
t > t(Z) means t — 1 > t(Z) and hence ^(X, Ox(Ft-i(Z))) = 0, so taking global sections 
of the exact sheaf sequence -> O x (Ft-i(Z)) -> O x (Ft(Z)) -> Eo ® O x (F t (Z)) -> 
0, we see H°(X,O x (F t (Z))) surjects onto H°(E ,0 Eo ® O x {F t (Z))). But since E is 
isomorphic to P 1 , we know that Oe <E> O x (F t (Z)) is isomorphic to C>E (t) and that 
H°(E ,Oe {1)) ® H (E ,0 Eo {t)) -> H°(E ,0 Eo {t + 1)) is surjective, and hence that 
H°(E , O E0 {l)) ® H°{E , O E0 ® O x {F t {Z))) -> ^°(E , + I)) is surjective. Taking 

global sections (denoted by r) of — > Ox - >■ Cx(-Eo) ~" ^ Cs (l) ~" > 0, tensoring by 
V = H°(X, O x (F t (Z))) and mapping by multiplication, one obtains the diagram 

o -> r x (Ox)®u -> r x (Ox(^ ))®u -> ^(0^(1)) ®v -> o 
o ^ r x (0 x (F t (z))) ^ r x (o x (F t+1 (z))) ^ r So (o E(1 (t + 1)) ^ o 

in which the leftmost vertical map is obviously an isomorphism and the rightmost vertical 
map as we saw is surjective, so the snake lemma gives an exact sequence in which the 
cokernels of the outer vertical maps are 0, hence the cokernel cok(ji{F t {Z))) of the middle 
vertical map also vanishes; i.e., v t {Z) = for t > t(Z) + 1.) 

In addition to Dubreil's bounds in Theorem II. 1.1, and Campanella's upper bounds 
v t +i(Z) < A 2 hz(t + 1) — e t mentioned in the proof above, Campanella also gave the lower 
bound that v t (Z) > ma,x{ A 3 h z (t), e' t }, with e'^ z ^ = 1 and e' t = otherwise (these bounds 

of course follow from v t {Z) — s t (Z) = A 3 /i^(t); Campanella was actually working in the 
more general situation of perfect codimension 2 subschemes of any projective space). 

Here is a result (a slight restatement of Lemma 4.1, [H7]) that in many cases turns out 
to be an improvement on the bounds above of Dubreil and Campanella, proved in a way 
very similar to the proof of Theorem II. 1.1 given above. It underlies many of the results 
of [H7], [H8], [HHF] and [FHH] . 

Lemma II. 1.2: Let Z = m\p\ H — • + m n p n be a fat points subscheme of P 2 with distinct 
points pi such that mi > 0, and let Z" = Z — p\ and Z' = Z +p\. Then max(hz(t + 1) — 
3h z (t) + hz»(t-l),0) < v t +i{Z) < h z (t + 1) - 3h z (t) + h z »(t - 1) + h z >{t). 
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II. 2. Exact Results 

Campanella's and Dubreil's bounds hold for any Z = m\p\ + ■ ■ ■ + m n p n C P 2 , not 
just when the points pi are general. Thus it is not surprising that the bounds are not 
always exact. For example, for Z = 3(p± H — ■ +p$) with pi general we have vg(Z) = 2 but 
max{A 3 /i z (8), ^} = 1 while A 2 h z (8) - e 7 = 3. 

Thus we can try to obtain exact results. The typical pattern for work on fat points has 
been first to obtain results when either the multiplicities rrii are small or the number n of 
points is small, and this is what we see regarding v. In particular, for Z = p\ + • • • +p n with 
Pi general, it is easy to see that oe(Z) is the least t such that t 2 + 3t+2 > 2n and that r(Z) is 
the least t such that t 2 + 3t+2 > 2n. And, as always, vt(Z) = unless oe(Z) < t < r(Z) + l, 
with, as always, v a (z)(Z) = h z (a(Z)), so for Z = p\ + ■ ■ ■ + p n only v a (z)+i(Z) remains 
to be found, and Geramita, Gregory and Roberts [GGR] proved for such subschemes of 
general points of P 2 with multiplicity 1 that [i a {z) has maximal rank (i.e., fjL a (z) is either 
surjective or injective, and hence v a (z)+i(Z) = rnax.{hz{oi{Z) + l) — 3hz(oi(Z)), 0}). Using 
different methods Ida [Id] has extended this to the case that Z = 2(pi + • • • + p n ) with pi 
general and n > 9. 

Results have also been obtained for large m if n is small. The first such I am aware of 
is that of Catalisano [Cat2], who determines vt(Z) f° r & h t an d any Z = m\p\ + ■ ■ ■ + m n p n 
as long as the points pi lie on a smooth plane conic; in particular, this handles the case of 
any Z involving n < 5 general points. 

To discuss Catalisano's result in more detail, let n be any positive integer (not nec- 
essarily 5 or less) and consider Z = m\p\ + ■ ■ ■ + m n p n for any n distinct points pi. If 
t > a(Z), let f t be a common factor of greatest degree for the elements of I(Z) t (i.e., 
ft = defines the fixed divisor of the linear system of curves given by I(Z) t ). Now let 
Z' t be ^(max{mj, mult Pi (ft)})Pi', thus f t spans I(Z' t )d, where d = Oi(Z' t ) is the degree 
of ft- In fact, d and Z' t can be found without finding ft- d = a(Z' t ) where among all 
Z" = m'{pi + • ■ ■ + m'nPn with < m" < m t and hz(t) = hz-z"(t — a(Z")) we choose Z[ 
to be that Z" for which J2i( m i ~ m i) ls least. 

In Catalisano's situation, the points pi are assumed to lie on a smooth conic. By 
[HI], hz and Z[ were already known and easy to compute for points on a conic, and, as 
noted above, it is enough to determine vt{Z) for t > a(Z). Catalisano's result, although 
expressed in her paper [Cat2] rather differently, can now be stated: 

Theorem II. 2.1: Let pi,...,p n be distinct points on a smooth plane conic, let Z = 
mipi + ■ ■ ■ + m n p n be a fat points subscheme of P 2 and let d = a(Z' t ), where Z[ is defined 
as above. Then for each t > a(Z) we have u t +i(Z) = hz(t + 1) — h z _ z ^(t + 1 — d). 

For a proof in a slightly more general case (the conic need not be smooth, for example, 
and the points can be infinitely near), see [H6]. As an aside, note that Theorem II.2.1 shows 
that u t+1 (Z) > only if a(Z) -l<t< /3(Z), since for t > f3(Z) we have f t = 1, so d = 
and Z' t = 0. 

Since n > 5 general points do not lie on a conic, Theorem II. 2.1 does not apply for 
n > 5 general points. Nonetheless, the inequality v t +i(Z) > h z (t + 1) — h z _ z > t (t + 1 — d) 
always holds (see Lemma 2.10(c), [H6]), although equality can fail for n > 5 general points 
since vt+i can be positive even if t > (3(Z). A complete solution for the case of any 
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Z = rriipi + • • • + m n p n for n general points pi was given for n = 6 by Fitchett [F2] , for 
n = 7 by me [H8] and finally for n = 8 by Fitchett, me and Holay [FHH]. 

Given Z = m\p\ + ■ ■ ■ + m n p n with pi general, the main result of [F2] is that n(F t (Z)) 
has maximal rank as long as F t (Z) is nef and n < 6. (A divisor D is nef if [D] ■ [H] > 
for every effective divisor H.) Given a divisor on the blow up X of P 2 at n general points, 
the main results of [H8] give an algorithmic reduction of the problem of determining the 
rank of ^(F) to the case that F is ample, and shows if n < 7 that (i(F) is surjective as 
long as F is ample, thereby solving the problem of resolving I(Z) for n < 7. (A divisor D 
is ample if [D] ■ [H] > for every effective divisor H.) 

However, if n = 8, reduction to the ample case is not enough, since examples show that 

fi(F t (Z)) can fail to have maximal rank even if F t (Z) is ample; take Z = 4(pH \-pr)+P8 

with t = 11, for instance (case (c.ii) of Theorem II. 2. 2). The main result of [FHH] boils 
down to giving a formula in nice cases together with an explicit algorithmic reduction to 
the nice cases. We now give a slightly simplified statement of the main result of [HHF]. 

Recall that an exceptional curve on a smooth projective surface S is a smooth curve 
C isomorphic to P 1 such that [C] 2 = —1 in Cl(S'). Assuming n = 8, let E x denote the set 
of classes of exceptional curves on X and for each exceptional curve C define quantities 
Ac and Ac as follows: For C = Ei for any z, let Ac = Ac = 0. Otherwise, let mc be the 
maximum of C ■ E\, . . . , C ■ E n , define Ac to be the maximum of mc and of (C ■ L) — mc 
and define Ac to be the minimum of mc and of (C • L) — mc- We then have: 

Theorem II. 2. 2: Let X be obtained by blowing up 8 general points of P 2 and let 
[Eq], [Ei], . . . , [Eg] be the associated basis of the divisor class group of X . Consider the 
class F = t[E ] — mi [Ei] — ■ ■ ■ — ms[E 8 ], where mi > • • • > m$. 

(a) If F ■ C > Ac for all C G "Ex, then n(F) has maximal rank. 

(b) If F C < Ac for some C e E x , then ker(p(F)) and ker(p(F — C)) have the same 
dimension. 

(c) If neither case (a) nor case (b) obtains, then either 

(i) F ■ {Eq — Ei — E 2 ) = 0, in which case cok(p, F ) has dimension h x (X, Ox(F — 
(E - Ei))) + h\X, O x (F - (E - E 2 ))), or 

(ii) [F] is [3E -Ex E 7 ] + r[8E - 3E ± 3E 7 - E 8 ] for some r > 1 

(in which case dim cok(/i(F)) = r and dim ker((i(F)) = r + 1), or 

(in) fi(F) has maximal rank. 

This theorem leads directly to an algorithm for computing resolutions for fat point 
subschemes Z involving n < 8 general points of P 2 . The MACAULAY 2 script f indres 
included at the end of this paper implements this algorithm to compute the values of 
vt(Z) and hz(t) for all a(Z) < t < r(Z) + 2. Since it does not rely on Grobner basis 
computations, it is in comparison quite fast. 

II. 3. The Quasi-uniform Resolution Conjecture 

What to expect for n > 8 remains mysterious. Whereas (as discussed below) by taking 
into account effects due to exceptional curves there results a reasonable conjecture for hz 
for any Z involving general points, doing the same for resolutions is harder. (For a partial 
result in this direction, see [Fl], which at least sharpens the bounds Campanella has given 
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on ft{Z). Also see Theorem 5.3 of [H7], which shows in a certain sense that behavior 
in the n > 8 case is simple asymptotically, and that it is the case of relatively uniform 
multiplicities that is not understood.) In fact, the results of [FHH] for n = 8 suggest that 
taking into account the exceptional curves may not be enough. Thus it is still unclear how 
vt(Z) should be expected to behave in general. 

If one puts a mild condition on the coefficients m^, however, there is reason to hope 
that behavior may be quite simple. In particular, say that Z = m\p\ + ■ ■ ■ + m n p n is 
quasi-uniform if: the points pi are general; n > 9 and mi = mg; and the coefficients 
mi > m2 > • • • > rn n > are nonincreasing. We then have the following Quasi-uniform 
Resolution Conjecture ([HHF]): 

Conjecture II. 3.1: If Z is quasi-uniform, then hz(t) = max{Pz(t),0} and Vt(Z) = 
max{h z (t) - 3h z {t - 1),0} (or, equivalently, v t (Z) = max{0, A 3 h z (t)}) for all t > 0. 

Assuming this conjecture, we can write down an explicit expression for the resolution 
of I(Z) for a quasi-uniform Z, as follows (see [HHF]): 

-> R[-a - 2} d © R[-a - l] c -> R[-a - if © R[-a] a -> I{Z) -> 0, 

where a = a(Z), a = hz(ot), b = max{/i^(« + 1) — 3hz(ct), 0}, c = max{— hz(ot + 1) + 
3hz(ot), 0}, and d = a + b — c — 1. 

Most of the evidence for this conjecture currently is for the uniform case (i.e., Z is 
quasi-uniform with all multiplicities mi equal to some single m) and mostly even then for 
small multiplicities. For example, the conjecture is true (and easy) for n = 9 (see [H6]), 
and whenever m < 2 (see [GGR] and [Id]). However, at the time that this is written, there 
are two situations in which the conjecture is known for unbounded multiplicities. The first 
is when n is a power of 4, and m > {\/n — 2)/4, in which case Conjecture II. 3.1 is true (in 
characteristic 0) by [HHF] and [Ev]. For the second (again in characteristic 0), see [HR], a 
typical example of which is n = d(d + 1) and m = d + 1 for any even integer d > 2. In this 
case, the bounds on a and r given by the modified unloading method due to Roe and me 
(using r = d(2d+ l)/2; see Section IV.1.3 and the end of Section IV.2) show a(Z) > t(Z) 
for Z = m(pi + ■ ■ ■ +p n ), from which it follows that I(Z) is generated in degree a(Z) and 
hence that the minimal free resolution of I(Z) is 

-> R[-d 2 -2d- l] d ( d + 2 ) -> R[-d 2 - 2d}^ 2 -> I(Z) -> 0, 

in agreement with Conjecture II.3.1. Additional evidence of various kinds for the conjecture 
is given in [HHF]. 

In order to facilitate checking this conjecture and exploring the problem of under- 
standing resolutions when Z need not be quasi-uniform, it is helpful to be able to compute 
resolutions directly. Although we are interested in general points, it is easiest instead to use 
random choices of points, with the expectation that this will usually give points that are 
general enough. It is possible to implement such a calculation very simply in MACAULAY. 
Here is an example of such a MACAULAY 2 script, provided to me by Hal Schenck, for 
computing resolutions of ideals I(^2 i miPi) for random choices of points pi G P 2 : 
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R=ZZ/31991 [x_0. .x_2] 

mixer = (l)->(i:=0; 

b:=ideal (matrix {{1}}**R) ; 
scan(#l, i->( 

f :=random(R~l,R~{-l}) ; 

g:=random(R~l ,R~{-1}) ; 

I:=(ideal (f I g))~(l#i); 

b=intersect (b , I) ) ) ; 
print betti res coker gens b; 
b) 

— Return the ideal of mixed multiplicity random fatpoints. Input is a list 

— with the multiplicities; e.g. mixer ({1 ,2,3}) returns the ideal of 

— I(pl)~l \cap I(p2)~2 \cap I(p3)~3, where pj is a (random) point in P~2, 

— and prints the betti numbers of the resolution 

— HKS 4/28 

By a slight modification, below, this script can be made to handle the uniform case 
(i.e., n random points each taken with the same multiplicity m): 

R=ZZ/31991[x_0. .x_2] 

unif = (n,m)->(« n « " points of multiplicity " « m « ":" « endl; 
i:=0; 

b:=ideal (matrix {{1}}**R) ; 
while i < n do ( 

f :=random(R~l,R~{-l}) ; 

g:=random(R~l,R~{-l}) ; 

I:=(ideal (f I g))~(m); 

b=intersect (b , I) ; 

i=i+l) ; 

print betti res coker gens b; 
b) 

— Example: unif ({3, 2}) returns the ideal of 

— I (pi) "2 \cap I(p2)~2 \cap I(p3)~2, where pj is a (random) point in P~2, 
— and prints the betti numbers of the resolution 
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III. Hilbert functions 

Nagata, in connection with his work on Hilbert's 14th problem, began an investigation 
of the Hilbert function hz for fat points subschemes Z = m\pi + - ■ ■ + m n p n with pi general, 
although his work was written up from the point of view of divisors on blow ups of P 2 (see 
[N2]). 

111.1. Nagata's Work 

In brief, Nagata in [N2] determines hz(t) (and thus a(Z)) for all t for any Z = 
m\p\ + ■ ■ ■ + m n p n with pi general as long as n < 9. In [Nl], he poses the following 
conjecture, which remains open unless n is a square, in which case Nagata verified it: 

Conjecture III. 1.1: Let Z = ui\p\ + ■ • ■ + m n p n for n > 9 general points pi G P 2 . Then 
a(Z) > (mi H \-m n )/y/n. 

Also implicit in [N2] is a lower bound (see (*), Section III. 2) for the values of the Hilbert 
function hz of I(Z). An easier lower bound comes from the fact, as discussed above, that 
h z (t) =P z {t) + h}{X,Ox{F t {Z))). Since h\X, O x (F t (Z))) > 0, it of course follows that 
hz(t) > max{P^(t), 0}. However, easy examples show that hz(t) > max{P^(t),0} can 
sometimes occur; in all such examples for which hz(t) is known, the difference hz(t) — 
max{P^(t), 0} has a geometric origin, being always precisely what one gets by taking into 
account exceptional curves. Taking the exceptional curves into account gives the more 
refined bound (*). 

To explain this, let X be obtained by blowing up n distinct points pi of P 2 . We 
have, as discussed above, the basis [Eq], . . . , [E n ] of the divisor class group Cl(X) of X. 
Because we are mostly interested in the case of n general points, technical issues force us 
to use the following definition. Let us say that an element v = Yli a i[Fi\ of C1(A) is an 
exceptional class if for general points pi there is an exceptional curve C C X with v = [C]. 
(The problem is that there may be no nonempty open set U of points (pi, . . . ,p n ) e (P 2 ) n 
for which all exceptional classes are simultaneously classes of exceptional curves, even 
though each exceptional class v is the class of an exceptional curve for some nonempty 
open U v C (P 2 ) n .) 

Nagata [N2] determined the set £{n) of exceptional classes. It turns out that 6(0) is 
empty, 5(1) = {[^i]}, and £(2) = {[E^, [E 2 ], [E -E 1 -E 2 ]}, while for n > 3 the set £(n) 
is the orbit W n [E n ] with respect to the action of the group W n of linear transformations 
on C\(X) generated by all permutations of {[-Ei], • • • , [E n ]} and (if n > 3) by the map 7 
for which 7 : [Ei] i-> [E t ] for % > 3, 7 : [E,] i-> [E ] - [E t ] - [E 2 ] - [E 3 ] + [E,] for < % < 3 
and 7 : [E ] 1— > 2[E ] — [Ei] — [E 2 ] — [E3]. (The map 7 can be regarded as a reflection 
corresponding in an appropriate sense to a quadratic Cremona transformation centered at 
Pi, p 2 and P3. The fact that W n is a reflection group was recognized by Du Val [DuV2], 
and extended and exploited by Looijenga [L].) 

111. 2. A Decomposition and Lower Bound 

For simplicity, assume n > 3. This will not be a serious restriction, since cases n < 3 
are easy to handle ad hoc, and in any case there are natural inclusions £(n) C £(n+l) for all 
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n, so a given value of n subsumes smaller values. Now Let \1/ be the subsemigroup of C\(X) 
generated by E(ri) and by the anticanonical class — Kx = [3Eq — E\ — ■ ■ - — E n ] of X. (With 
respect to the action of W n on C\(X), \1/ is essentially Tits' cone [Ka]; thus there exists a 
fundamental domain for the action of W n on For any F G ^, it turns out that there is 
a unique decomposition F = Hp + Np with (dropping the subscripts) H,N G ^ such that 

H • v > for every exceptional class v, H ■ N = 0, and either N = or N = a\V\-\ \-a r v r 

for some exceptional classes vi and integers a\ > 0, such that Vi ■ Vj = for all i ^ j. (It 
is easy to compute this decomposition. By recursively applying 7 and permutations in a 
straightforward way, for any F G Cl(X) one can find an element w G W n such that either 
wF ■ [E ] < 0, or wF ■ [E - E{\ < 0, or such that wF = a [E ] + J2i>o a-i[Ei] with a > 0, 
flo + a i + a 2 + ^3 > and a\ < • • • < a n . But if either wF ■ [E ] < or wF ■ [E — E{\ < 0, 
then F^ty, while otherwise there are two cases. Either > wF ■ [Eq — E\ — E2] = 
ciq + ai + a,2, in which case H = w _1 [(2ao + a± + 02)^0 — (^o + a>2)Ei — (ao + ai)^] and 
N = w- 1 ((-ai - a 2 - ao) [E - E 1 - E 2 ] + E a , >o a * ) » or wF ' [^0 - £?i - £? 2 ] > and we 
have H = w~ 1 (ao[E Q ] + Yla <o a i[Ei\) and = w~ 1 (J2 a >0 ai[E i \). An implementation 
of this procedure is given by the script decomp provided in this paper.) 

It is true (and more or less apparent from [N2]) for general points pi that if 
h?{X,O x {F)) > then F G hence F = H + N as above, and h°(X, O x (F)) = 
h°(X, O x {H)) > (H 2 — K ■ H)/2 + 1. For any F G C\(X), define e(F) to be unless 
F G in which case set e(F) to be the maximum of 1 + (Hp — K ■ H F )/2 and 0. We then 
get the lower bound 

h°(X,O x (F))>e(F). (*) 

(The script homcompdim included at the end of this paper computes e(F). For F = 
d[E ] — mi [-Ei] — ••• —m n [E n ], we have e(F) = homcompdim({<i, {mi, m n }}), which, if 
the multiplicities m^ are nonnegative, is also thus a lower bound for the dimension of the 
homogeneous component of I(Z) of degree d for Z = mipi + ■ • • + m n p n .) 

III.3. The SHGH Conjectures 

It follows from Nagata's work (see Theorem 9, [N2]) that in fact h°(X, O x {F)) = e(F) 
for n < 9 general points. What occurs for n > 9 is not known, but I [HI] (also see [H3]), 
Gimigliano [Gil] (also see [Gi2]) and Hirschowitz [Hi] independently gave conjectures for 
explicitly computing h°(X, Ox{F)) for any n. These conjectures are all equivalent to the 
following conjecture, which states that e(F) is the expected value of h°(X, Ox(F)): 

Conjecture III. 3.1: Let X be the blow up of n general points of P 2 and let F G Cl(X). 
Thenh°(X,O x (F))=e(F). 

It is interesting to compare Conjecture III. 3.1 with an earlier conjecture posed by Segre 
[Seg], giving a conjectural characterization of those classes F such that h°(X, O x (F)) > 
max{0, (F 2 - K ■ F)/2 + l}: 

Conjecture III. 3. 2: Let X be the blow up of n general points of P 2 and let F G C1(X). 
Ifh°(X, O x {F)) > max{0, (F 2 — K ■ F)/2 + 1}, then the fixed locus of \F\ has a double 
component. 
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It is easy to show that Conjecture III. 3.1 implies Conjecture III. 3. 2; the fact that 
Conjecture III.3.2 implies Conjecture III. 3.1 is essentially Theorem 8 of [N2]. Thus I will 
refer to these conjectures (in any of their forms) as the SHGH Conjecture. Since Nagata's 
paper is hard to read, the equivalence of Conjecture III. 3.1 and Conjecture III. 3. 2 was only 
recently recognized and proved by Ciliberto and Miranda. Here is a sketch of a proof. 

Theorem III. 3. 3: Conjecture III. 3.1 is equivalent to Conjecture III. 3. 2. 

Sketch of proof: To see that Conjecture III. 3.1 implies Conjecture III. 3. 2, assume that 
h°(X,O x (F)) > max{0, (F 2 - K ■ F)/2 + 1} for some F. Thus h°(X,O x (F)) > 0, 
so F £ f and hence we have a decomposition F = H + N, as described above, with 
N = a±vi + ■ • ■ + a r v r for some exceptional classes Vi and ai > 0. By Conjecture III. 3.1, 
we have < h°(X, O x {F)) = e(F) = (H 2 - K ■ H)/2 + 1, and by substituting H + N in 
for F we see that (F 2 - K ■ F)/2 + 1 = (H 2 - K ■ H)/2 + 1 unless a* > 1 for some i, in 
which case Vi is the class of a curve occurring (at least) doubly in the base locus of \F\, 
proving Conjecture III. 3. 2. 

Conversely, assume Conjecture III. 3. 2. Among all F for which h°(X, Ox(F)) — e(F) 
fails, choose one having as few fixed components as possible (i.e., for which the sum of 
the multiplicities of the fixed components is minimal). As before we have F = H + N, 
but N = by minimality (since h°(X,O x (H)) = h°(X,O x (F)) > e(F) = e{H)). Since 
F = H, by construction of H we have F ■ E > for every exceptional class E. 

Now say some reduced irreducible curve C occurs as a fixed component of \F\ with 
multiplicity at least 2. Thus h°(X,O x (2C)) = 1, hence h°(X,O x (C)) = 1 and by 
Conjecture III.3.2 we have (C 2 - C ■ K)/2 + 1 = 1 so C 2 = C ■ K. Therefore the 
genus g c of C is (C 2 + C ■ K)/2 + 1 = C 2 + 1; i.e., C 2 > -1. On the other hand 
1 = h°(X, O x {2C)) > {AC 2 - 2C ■ K)/2 + 1 = C 2 + 1 so C 2 < 0. 

If C 2 = —1, then gc = 0, so C is an exceptional curve. From — > O x (F — C) — > 
O x (F) -> O c (C ■ F) -> it follows that h 1 (X,O x (F)) = (since C ■ F > implies 
h x (X,O c {C ■ F)) = 0, while h 1 (X,O x (F - C)) = by minimality), which contradicts 
failure of h°(X, O x {F)) = e(F). 

If C 2 = 0, then g c = 1, so C is an elliptic curve. From -> O x (F - C) -> Cx(^) -> 
Cc <8> ^(i 71 ) -> it follows that ^(X, Cx^)) = (as long as we see h x {X,O c ® 
Cx(-f 1 )) = 0, since /i 1 (X, O x (F — C)) = by minimality), which again contradicts failure 
of h (X, O x (F)) = e(F). But C (being irreducible of nonnegative selfintersection) is 
obviously nef, so C-F > 0. Since C is elliptic, C F > guarantees h x {X, O c ®O x (F)) = 
unless the restriction Oc <8> O x {F)) of Ox(-f) to C is trivial. But because the points 
blown up to obtain X are general, Oc <S> O x (F) cannot be trivial. (In fact, up to Cremona 
transformations, which is to say up to the action of W n , we can assume that [C] = [SEq] — 

[Ei H h Eg] . For Oc <S> O x (F) to be trivial we would need F ■ C = and, writing [F] as 

[dEg — m\E\ — • m n E n ], we would have 3d — (mi H — ■ + mg) = 0; i.e., [dEg — m\E\ — 

mgEg] is the class of an effective divisor perpendicular to [3Eg] — [E\ H h-Eg], but 

for general points pi, . . . ,pg the only such classes are multiples of [3Eg] — [E\ + ■ ■ • + Eg] 
itself, for which it is easy to check the restrictions to C are not, in general, trivial.) {> 
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III. 4. Evidence 

It is worth mentioning that it is not hard to show (see [HHF]) that the SHGH Conjec- 
ture implies that hzit) = max{P^(t), 0} if Z is quasi-uniform, which is part of Conjecture 
II. 3.1 posed above. In particular, if F = d[Eo] — m[Ei + ■ ■ ■ + E n ] where the E± are ob- 
tained by blowing up n > 9 general points of P 2 , then Conjecture III. 3.1 predicts that 
h (X, Ox(F)) equals the maximum of and (F 2 — K ■ F)/2 + 1. Proving this equality is 
trivial if m = 1, and was proved for m = 2 by Alexander and Hirschowitz in a series of 
papers culminating in [AH1] (it is worth noting that these papers address P^ for all N). 
More generally, given any positive integer M and M > mi > for all i, [AH2] shows for 
any Z = m\p\ + ■ — h m n p n (in any projective space) that hzit) = max{Pz(t), 0} for all t 
as long as n is sufficiently large compared with M. More explicitly, Ciliberto and Miranda 
[CM1], [CM2] have verified the SHGH Conjecture in characteristic for all m < 12 for 
any n > 9 (see also [Sei]), and Mignon [Mi] has now verified the SHGH Conjecture for all 
n > 9 for any F = d[Eo] — m\ \E\\ — • • • — m n [E n ] as long as rrii < 4 for all i. 

Whereas all of the explicit verifications of the SHGH conjecture described above as- 
sume multiplicities at most 12, two methods now exist that work for multiplicities which 
in some cases can be arbitrarily large; both assume that the characteristic is 0. The first 
is the recent result of Evain [Ev], which, for example, shows that h°(X, Ox{F)) equals the 
maximum of and (F 2 — K ■ F)/2 + 1 for any F = d[Eo] — m\E\ + • • ■ + E n ] as long as 
X is obtained by blowing up n general points with n being a power of 4. The second is 
the modified unloading method [HR] jointly due to me and J. Roe (see Section IV. 1.3 and 
the end of Section IV. 2), which gives very tight bounds on a and r. With good enough 
bounds, one can sometimes show a(Z) > r(Z), but anytime one knows a(Z) > t(Z) it 
immediately follows that the SHGH Conjecture holds for Z. In fact, there are numerous 
examples for which the bounds from [HR] are good enough to show ol(Z) > t(Z) and 
hence that the SHGH Conjecture holds, including certain infinite families of examples 
Z = m(pi + ■ ■ ■ + p n ) such as with n = d(d + 1), m = d + 1 for any even integer d > 2 
(mentioned in Section II. 3) or with n = d 2 + 2 and m = d(d 2 + 1) + d(d + l)/2 for any 
d > 2 (see Corollary V.2 of [HR] for these and other examples). 

We close this section with the comment that the script f indhilb computes the SHGH 
conjectural values of the Hilbert function of I(mipi + • • • + m n p n ) for general points pi. 
For n < 10, these values are the actual values. No separate script for the case of uniform 
multiplicities is included since for Z = m(pi + ■ ■ ■ + p n ) for n > 9 general points pi, the 
conjecture is simply that hzit) is the maximum of and (t 2 + 3t + 2 — nm{m + l))/2. 

IV. Bounds 

Rather than trying to prove the SHGH Conjecture directly, a good deal of work has 
been directed toward obtaining better bounds on a and r. The values of a and r predicted 
by the SHGH Conjecture give upper and lower bounds, respectively; in particular, a(Z) is 
less than or equal to the least t such that e(F t (Z)) > 0, and r(Z) is greater than or equal 
to the least t > such that hz(t) = Pz(t)- Thus what is of most interest are lower bounds 
on a and upper bounds on r. Bounds on a are especially of interest, since a sufficiently 
good lower bound on a may equal the upper bound (and presumed actual value) of a given 
by the SHGH Conjecture, and, as discussed above, if a always has its conjectured value 
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then the full SHGH Conjecture is true. 

Unfortunately, such tight bounds are so far fairly rare, but there are some, such as 
n = d(d + 1) points taken with multiplicity m = d + 1 as discussed above, for which (by 
precisely this method of tight bounds) the Hilbert function and resolution are known. For 
two additional examples, consider F = d[E ] — m[Ei + ■ ■ ■ + E n ] where X is obtained by 
blowing up n general points of P 2 with n being either 16 or 25. Although Evain's method 
handles these cases, at least in characteristic 0, an alternate approach is to notice that the 
inequality h°(X, O x (F)) > (F 2 -K-F)/2 + l with F = d[E ] -m[E 1 + - • - + E n ] guarantees 
that a < m\^fn \ + |~( \\fn \ — 3)/2] . For n = 16 general points this gives a < 4m + 1 while 
for n = 25 this gives a < 5m + 1. But Nagata's result [Nl] that a > rriy/n when n is a 
square bigger than 9 now shows that a = 4m + 1 for n = 16 and a = 5m + 1 for n = 25. 
By [HHF], r = a in these cases, which determines hz for all t (and even the resolution of 
I(Z) when n = 16). 

Some of the bounds discussed below are algorithmic in nature, and hard to give simple 
explicit formulas or estimates for. Thus, to compute them, I have included at the end of 
this paper two MACAULAY 2 scripts, bounds (/) and unif bounds (/); in the former case 
I = {mi, . . . , m n } (corresponding to taking n general points with multiplicities mi, . . . , m n ) 
while in the latter case / = {n, m} (corresponding to taking n general points each with 
multiplicity m). 

IV. 1. Bounds on a 

By Nagata's work [N2], the exact value of a(Z) is known for any Z = mipi + - ■ -+m n p n 
with pi general and n < 9, and in such cases can be computed by running the script 
f indalpha or unif f indalpha. For Z = m(pi + ■ ■ ■ + p n ), n < 9, it is easy to be explicit: 

a(Z) = \c n m\, where ci = C2 = 1, C3 = 3/2, 04 = 05 = 2, cq = 12/5, cj = 21/8, c§ = 48/17 
and eg = 3. 

For n > 9, f indalpha or unif f indalpha only give upper bounds for a, although 
the upper bounds given should be, according to the SHGH Conjecture, the actual values. 
Thus most interest is in finding lower bounds on a, and a number of such have been given. 
Let n > n' > 9 and consider Z = m(pi + ■ ■ ■ + p n ) and Z' = m(pi + • • • + p n '), where 
the points pi are general. It is easy to see that a(Z) > a(Z'). Since Nagata [Nl] proves 
that a(Z') > m\fnJ if n' is a square, it follows (taking n' = [^/n] 2 when n is 16 or more) 
that a(Z) > mLv^J- A complete proof is somewhat tricky; we treat the slightly weaker 
inequality a(Z) > m[y/n\ in the next section. 

IV. 1.1. Bounds by testing against nef divisors 

The inequality cx(Z) > m\_^/n\ follows easily (for any n) by specializing [v^J 2 °f the 
points pi to a smooth plane curve C of degree [v^-J • The class C of the proper transform 

of C to the blow up X of P 2 at the points pi is [v^J [-^o] — {[E\ + V E^^ 2 ]), which 

is nef, but a(Z)[E ] — m[E\ + h E n ] is (by definition of ot(Z)) the class of an effective 

divisor, so the intersection C ■ (a(Z)[Eo\ — m\E\ + ■ ■ ■ + E n ]) = a{Z)\_^Jn\ — m[y/n\ 2 is 
nonnegative, which gives oe(Z) > mY^/n\. More generally, the same argument works for 
Z = mxpx H h m n p n , giving a(Z) > (mi H h m^p)/ [^/n\ . 

Alternatively, by specializing all n points to a curve of degree \y/n\ , the same argument 
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(using the fact that now \y/n \ [E ] — [E\ + ■ ■ ■ + E n ] is nef) gives the inequality a(Z) > 
mn/\y/n) for Z = m(p\ + ■ ■ ■ +p n ), and a(Z) > (mi + • • ■ + m n )/[^fn\ for Z = m\p\ + 
• • • + m n p n . More generally, we have the following extension of the main result of [H9]: 

Theorem IV. 1.1.1: Let Z = ui\p\ + ■ ■ ■ + m n p n for general points pi G P 2 with n > 1 
and mi > • • > Tn n , and let r < n and d be positive integers. Given nonnegative 
rational numbers (not all 0) ciq > a± > ■ ■ ■ > a n > such that a^d 2 > a\ + ■ ■ ■ + a r and 
ra > ai H h a n , then a(Z) > a i m i) / ( a od) ■ 

Sketch of proof: Note that by multiplying by a common denominator, we may assume 

that each aj is a nonnegative integer. Consider the class F = [aodE — a\E\ a n E n ] 

on the surface X obtained by blowing up the points pi. First, specialize (as in the proof of 
the main result of [H9]) to certain infinitely near points; in particular, such that [Ei — Ei + i] 
for each < % < n is the class of an effective, irreducible divisor on the specialization X' 
of X , and such that d[Eo] — \E\ + ■ ■ ■ + E r ] is the class C of the proper transform of a 
smooth plane curve. Now F is nef on X' and hence on X. To see this, note that: F ■ C > 
since ao<i 2 > a\ + ■ ■ ■ + a r ; F ■ [Ei — E i+ i) > for alH > since aj > aj + i; F ■ E n > 
since a n > 0; and F is a nonnegative integer sum of the classes C, [Ei — i^+i], i > 0, 
and [E n ] since ao > a\ and rao > a\ + • • • + a n . Thus F is a sum of effective classes (in 
particular, of ciqC and various multiples of the [Ei — Ei+i] and E n ), each of which it meets 

nonnegatively; thus F is nef and so F meets a(Z)E — m\E\ m n E n nonnegatively, 

from which our result follows. <0> 

Finding an optimal bound for a given Z using Theorem IV. 1.1.1 involves solving a 
linear programming problem (note that we may normalize so that ao = 1), not to mention 
the problem of identifying the best choices of r and d. In case the multiplicities are 
all equal, it is not hard to show that optimal solutions (for given r and d) to this linear 
programming problem are given in parts (a) and (b) of the following corollary. These need 
not always be optimal if the coefficients are not all equal, so we consider in parts (c) and 
(d) some additional possibilities. 

Corollary IV. 1.1. 2: Let Z = m\p\ + h m n p n for general points pi G P 2 with n > 1 

and mi > ■ ■ ■ > m n > 0, let r < n and d be positive integers and let m be the mean of 
m 1 ,---,m n . 

(a) If r 2 > nd 2 , then ot(Z) > mnd/r. 

(b) Ifr 2 < nd 2 , then ol[Z) > mr/d. 

(c) Ifd 2 > r, then a(Z) >{m 1 ^ + m r )/d. 

(d) Assume d 2 < r and let j be an integer, < j < d 2 . 

(i) If j = 0, then a(Z) > (mi + • ■ ■ + m^)jd. 

(ii) If j > 0, let t = min{r + (r — d 2 )(r — d 2 + j) /j, n} and set m t +i = if t = n; 
then 

Sketch of proof: Each part of the corollary applies Theorem IV. 1.1.1 for various values 
of the ai. For (a), take ao = r and a\ = ■ ■ ■ = a n = d 2 . For (b), take ao = n and ai = r, 
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% > 0. For (c), take ai = 1 for % < r and ai = for i > r. For (d)(i), take = 1 for 
i < d 2 and = for i > d 2 . For (d)(ii), take a« = 1 for i < d 2 — j and = j/(r — <i 2 + j) 
for d 2 — j < i < [t\. If £ = n, then m t+ i = (and so is t — but if £ < n, then take 
a t+ i = (t- [t\)j/(r-d 2 + l). 

One can formally verify that the values of the given in (d)(ii) satisfy the necessary 
conditions to apply Theorem IV. 1.1.1, but it may be helpful to briefly discuss how these 
values come about. The idea giving rise to the values of in (d) (ii) is to find extremal 
sets (one set for each j) of values of the a^, with the hope that for any given Z one set will 
be close to an optimal solution that might be found by linear programming. By setting an 
equal to 1 (a normalization we clearly can always do), we bound the values of the other ai 
above by 1. Since the multiplicities rrij are nonincreasing, any optimal solution for the 
must also be nonincreasing. Intuitively, we would want to keep as many of the a« equal to 
1 as possible. But in order to satisfy d 2 > a± + ■ ■ ■ + a r we can keep at most the first d 2 
of the di equal to 1, in which case all of the other ai would have to be 0. Depending on 
the values of the m^, however, we may be better off if we can make enough of the other 
di positive. So, given j, we leave ai, . . . , a^-j alone, and spread a^-j+i, • • • , a^, which 
are each 1 to start with, evenly over a^2_ J+1 to a r , which reduces a^-j+i, . . . , a^ 2 from 

1 to j / (r — d 2 + j), and raises a^+i, . . . , a r from to j / (r — d 2 + j), while keeping the 
condition d 2 > a\ + ■ ■ ■ + a r satisfied at equality. Now, although this may have worsened 
things (since we may well have reduced a\m\ + ■ • ■ + a r m r ), we can hope to more than 
make up for this since we can now increase some of the remaining a« from (which they 
were before) to j / (r — d 2 How many of the ai which we can increase is limited by the 
condition r > a± + ■ — h a n = d 2 + a r +i + • • • + a n ; moreover, because of fractional effects, 
the last ai which we can manage to increase from might be limited to being increased 
only by a fraction of j / (r — d 2 + j), which accounts for the anomalous behavior of a t +i- <0 

The bounds given in Corollary IV. 1.1. 2 can be computed by running the scripts unif- 
bounds or bounds. The script ezbhalphaD, which is called by bounds, checks all possible 
r, d and j from Corollary IV. 1.1. 2(d). 

IV. 1.2. Bounds by unloading 

As an alternative to Theorem IV. 1.1.1, we can use a process that can conveniently be 
referred to as unloading. The idea is based on the fact that given a divisor class D on a 
surface X and some finite set S of classes of effective, irreducible divisors, if for some F G S 
we have F ■ D < 0, then clearly D is the class of an effective divisor if and only if D — F is. 
Unloading (in a sense that is slightly more general than its use in the literature) consists 
of checking D ■ F for each F G S, and replacing D by D — F whenever D ■ F < and 
continuing with the new D. (For the classical notion of unloading, see pp. 425-438 of vol. 

2 of [EC], where it is referred to as scaricamento, or see [DuVl].) Eventually, D reduces 
to a class D' such that either D' is obviously not effective (because, perhaps, D' ■ E < 
or D' ■ (E - E x ) < 0) or such that D' ■ F > for all F G S. 

With respect to the specialization used in the proof of Theorem IV. 1.1.1, we can take 
S to consist of the classes [Ei — E i+ i] for < i < n, [E n ] and d[E ] — [E\ + ■ • ■ + E r ], and 
we look for the largest t such that D = [tEo — (mi Si + • ■ • + m n E n )] unloads to a class D' 
with D' • (Eq — Ei) < 0, in which case t + 1 is a lower bound for a(Z). I have included the 
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script bhalpha to compute the bound obtained via unloading with respect to any chosen 
r and d. 

In the special case that r < d 2 , then d[Eo] — [E\ + ■ ■ ■ + E r ] is nef. It is not hard to 
then see that the result of the unloading process is the same as just testing against this nef 
divisor, hence, assuming that mi > • • • > m n , we get the bound a(mipi + ■ ■ ■ + m n p n ) > 
(mi + • • ■ + m r )/d. One can also give a formula for the result of this unloading process 
in another extremal case, 2r > n + d 2 . In this case, for Z = m(pi + • • • + p n ), we have 
a(Z) > 1 + ud + mm{d — 1, \p/d] — 1}, where u > and p are defined by mn = ur + p, 
with < p < r. 

The idea of bounding a using unloadings and specializations to infinitely near points 
is due to Roe [Rl], who actually uses a sequence of increasingly special specializations of 
infinitely near points, applying unloading after each specialization. Roe uses a sequence of 
n — 2 specializations, corresponding to sets Si, 3 < i < n, of classes of reduced irreducible 

divisors, where Si = {[E n ]} U {[Ej - E j+1 ] : 1 < j < n} U {[Ex - E 2 Ei]}. Starting 

with F t = tE — rriiEi — m 2 -E , 2 — ■ • • — m n E n , Roe's algorithm consists of unloading F t 
with respect to S3 to get F^ 3 \ then unloading F^ with respect to £4 to get F^ 4 \ etc., 
eventually ending up with F^ = tE — Ex — E 2 — ■ ■ ■ — m^E n . Roe's bound 
is then a(Z) > m["' , which comes from the fact that [F^] and hence [F t ] cannot be 

classes of effective divisors unless t > m^. This bound can be computed with the scripts 
unif roealpha and roealpha. 

Although it is hard to give a simple formula for the exact value of the result of this 
method, an asymptotic analysis by Roe [Rl] shows that his unloading procedure gives a 
lower bound for a(Z) which is always better than m{^n — 1— 7r/8), for Z = m(pi + - ■ -+p n ) 
with n > 2 general points. It should be noted however that this formula often substantially 
understates the result of the full algorithm. 

IV. 1.3. Bounds by a modified unloading 

Assume a specialization as in the second paragraph of Section IV. 1.2; we may assume 
d[Eo] — [Ei + • ■ • + E r ] is the class of a smooth curve C. Unloading F = [tEo — m\E\ — 

m n E n ] with respect to the set S of Section IV. 1.2 uses the fact that, if C ■ F < 0, then 

F — [C] is the class of an effective divisor if and only if F is. However, the requirement 
C ■ F < can be relaxed, since all we really need is h°(C, Oc{F)) = in order to ensure 
that F — [C] is the class of an effective divisor if and only if F is. By joint work with J. 
Roe [HR] , using the notion of a flex of a linear series on C, one can show (in characteristic 

0) that h°(C, O c {F)) = if either t < d and (t + l)(t + 2)/2 < mi H + m r , or t > d-3 

and F ■ C < (d — l)(d — 2)/2 — 1. (Recall that C is the proper transform of a plane curve 
C . The idea is to choose p± G C so that it is not a flex for the complete linear series 
associated to the restrictions to C of the divisors occurring during the unloading process. 
This is automatic in characteristic as long as p\ is a general point of C", but in positive 
characteristics every point of C may be a flex for a given, even complete, linear series [Ho].) 
Using this test in place of the more stringent test C ■ F < discussed at the beginning 
of Section IV. 1.2 gives what may be called the modified unloading procedure. Since this 
modified procedure uses a less stringent test, a larger (or at least as large) degree is needed 
to pass the test, so it gives bounds on a which are at least as good as the original unloading 
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procedure. These new bounds can be computed by running HRalpha or unif HRalpha. 

Although this modified unloading procedure is somewhat difficult to analyze in 
general, in two extremal cases Roe and I can derive the following simple bounds for 
Z = m(pi + • • • + p n ), for which define u > 0, p and s by requiring ran = ur + p, 
with < p < r, where s is the largest integer such that (s + l)(s + 2) < 2p: 

• if 2n > 2r > n + d 2 , then a(Z) > s + ud + 1. (The bound given by this formula 
can be computed by running ezunif HRalpha.) 

• if d(d+ 1)/2 < r < min{n, d 2 }, then a(n; m) > 1 + min{ \_(mr + g — l)/d\, s + ud}. 
(The bound given by this formula can be computed by running ezunif HRalphaB.) 

IV. 1.4. Bounds using \1/ 

The subsemigroup ^ C C\(X), introduced in Section III. 2, contains the subsemigroup 
of classes of effective divisors. Thus, given Z : the least t such that F t (Z) G ^ is a 
lower bound for a(Z). This sometimes gives an optimal bound. For example, if Z = 
90pi + 80p 2 + 70p 3 + 60p 4 + 50p 5 + 40(p 6 +p 7 +p 8 ) + 30p g + 20pi + 10pn, then the least t 
such that F t (Z) is in \1/ is 179, hence in fact a(Z) = 179, since e(Fi 79 (Z)) > 0. Finding the 
least t such that F t (Z) E \1/ is somewhat tedious, so I have provided the script Psibound 
for doing so. 

IV. 1.5. Comparisons 

For subschemes Z whose multiplicities are not too uniform, the lower bound on a(Z) 
given by testing against \1/ can be the best, as it is for Z = 90pi + 80^2 + 70j>3 + 60j>4 + 
50p 5 + 40(p 6 + pi + ps) + 30p 9 + 20pi + 10f»ii (see Section IV. 1.4). For example, Roe's 
method [Rl] of unloading gives a(Z) > 162, and the best result achievable using Corollary 
IV. 1.1. 2 turns out to be a(Z) > 173, whereas testing against \1/ gives a(Z) > 179 (and 
hence a(Z) = 179 as discussed above). 

However, if the multiplicities are fairly uniform, testing against \1/ does not give a very 
good bound. For example, for Z = m(pi + ■ ■ ■ + p n ) with n > 9, it is easy to see that 
F t (Z) e \& for all t > 3m, so testing against \P gives the bound a(Z) > 3m. This compares 
poorly with bounds via the other methods, which are typically very close to, but usually 
less than, m^fn. (Currently only the unloading method of [Rl] and the modified unloading 
method, discussed in Section IV. 1.2 and Section IV. 1.3, resp., ever are substantially better 
than rriy^n, and even these only when m is not too large compared to n.) 

Thus for uniform subschemes Z = m(pi + ■ ■ ■ + p n ) one is better off using some 
method other than testing against such as testing against nef divisors, as discussed in 
Section IV. 1.1. In this case one has, for any d and r, easy to implement tests, as given 
in Corollary IV. 1.1. 2. By comparison, the result of unloading with respect to the divisor 
C = dE — (Ei + ■ ■ ■ + E r ), as discussed in Section IV. 1.2, is, except in certain special 
cases, harder to compute since there is not always a simple formula for the result. Since 
one rarely gets something for free, it is not surprising, for given r and <i, that the bounds 
given by testing against a nef divisor are never better than those given by unloading. 

To see this, let Z = m(pi + ■ ■ ■ + p n ), and assume a(Z) > t is the bound given by 
unloading with respect to C. Also, with respect to the same r and d, let F = ao-Eo — 
(ai-Ei + • — h a n E n ) be the nef test class in the proof of Theorem IV. 1.1.1. The unloading 
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method unloads a divisor D = tE — m(E\ + ■ ■ ■ + E n ) to a divisor D' which meets C, 
E n and Ei — E^ + i, for all i, nonnegatively. But D — D' is a sum of multiples of these 
same divisors, which are all (linearly equivalent to) effective divisors, so each meets F 
nonnegatively. In addition, F is a sum of these same divisors, each of which D' meets 
nonnegatively, so D' ■ F > too. Thus D ■ F > 0, which shows that testing against the nef 
divisor can never rule out the candidate obtained by unloading. 

Moreover, if r 2 > nd 2 , unloading can definitely be better. For example, take n = 22 
and m = 3. Then the best choice of r and d with r 2 > nd 2 is r = 19 and d = 4, while the 
best choice of r and d with r 2 < nd 2 is r = 14 and d = 3. Using Corollary IV.1.1.2(a,b) 
with either choice of r and d gives a > 14, but unloading with respect to r = 19 and d = 4 
gives a > 15. Since Corollary IV. 1.1. 2 is optimal in this case, we see unloading sometimes 
gives a better result than can be obtained by any use of Theorem IV. 1.1.1. 

On the other hand, for Z = m(pi + • ■ ■ + p n ) with r 2 < nd 2 and r < n, the bound 
a(Z) > mr/d obtained by testing against a nef divisor, cannot be improved by unloading 

with respect to dE — {Ei H \-E r ), and hence unloading and testing against a nef divisor 

give the same result in these circumstances. (This is because for unloading to give a better 
bound, the class D = [\mr /d~\Eo — m(Ei + ■ ■ ■ + E n )] would have to unload to something 
obviously not effective, but unloading cannot get started unless D meets dEo — -Ei — • • ■ — E r 
negatively, which it does not.) But as the example of the preceding paragraph shows, if 
r 2 < nd 2 , although one cannot do better than mr/d by unloading with respect to r and 
d, one can still hope to do better than mr/d by unloading using some choices r' and d! in 
place of r and d. 

Since the modified unloading procedure of Section IV. 1.3 uses a less stringent test 
than does unloading, as in Section IV. 1.2, with respect to C = dE — (Ei + ■ ■ ■ + E r ) (in 
the sense that in order to be allowed to subtract C and continue the unloading process, 
for the former the intersection of with C can in most cases be as much as g — 1, where g 
is the genus of C, whereas for the latter the intersection must be negative), we see that 
bounds obtained via the latter method can never be better than those obtained by the 
former. The advantage of the latter method is that no hypotheses are required on the 
characteristic. 

There is also Roe's unloading method [Rl], discussed in Section IV. 1.2. As shown in 
[H9], for m sufficiently large compared to n, the results of Corollary IV. 1.1. 2 are always 
better than Roe's unloading method. However, when m is not too large compared with n, 
examples indicate that Roe's method gives the best bounds currently known. Consider, 
for instance, two examples using modified unloading (Section IV. 1.3). For n = 1000 and 
m = 13, Roe's method gives a > 421, whereas modified unloading, using r = 981 and 
d = 31, gives a > 424, and the SHGH conjectural value of a is 426. For n = 9000 and 
m = 13, things become reversed: Roe's method gives a > 1274, while modified unloading 
using r = 8918 and d = 94 gives only a > 1267; the SHGH conjectural value of a in this 
case is 1279. 

An interesting feature of these examples is that in both cases the bounds are better 
than \m^Jn\ + 1, conjectured by Nagata (Conjecture III. 1.1): For m = 13 and n = 1000, 
we have [wa/tTJ + 1 = 412, while for m = 13 and n = 9000, we have [/fty^J + 1 = 1234. 
Indeed, whereas most known lower bounds for a(Z) for Z = m(pi + ■ — Vp n ) are less than 
m\fn (since after all Nagata's conjecture is still open), the method of [Rl] and that of 
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modified unloading are among the few that in certain situations gives bounds that can be 
substantially better than rriy/n. In particular, if m is no bigger than about y/n, the method 
of [Rl] consistently (and probably always, although this looks hard to prove) gives a lower 
bound that is at least as big as my/n, and gets better as m decreases until, for m = 1 
it is easy to show that it gives the actual value of a(Z). If one chooses r and d carefully 
(depending on n), examples indicate that the modified unloading procedure does nearly as 
well as the method of [Rl] when m is small compared to n, and is substantially better for 
larger m. The method of [Rl], of course, has the advantage of being characteristic free and 
does not depend on careful choices of other parameters. The modified unloading method, 
on the other hand, sometimes gives a lower bound which is equal to the SHGH conjectural 
value (which is known to be an upper bound), and thus determines a exactly (as happens, 
for example, when n = d(d + 1) and m = d+lford>2 even, as discussed in Section II. 3, 
or n = 38 with m = 200, as mentioned in Section III. 4). 

Thus, in terms of getting the best bound for a given Z ', the modified unloading method 
(at least in characteristic 0) is often the best. It has, compared with methods (such as 
Corollary IV. 1.1. 2) which test against nef divisors, the disadvantage of being harder to 
compute, unless special values for r and d are chosen for which a formula applies. But 
since Corollary IV. 1.1. 2 works for essentially any r and d, sometimes one can do better by 
applying Corollary IV. 1.1. 2 than one can by applying the formula of Section IV. 1.3 where 
one's choices of r and d are more restricted. 

This raises the question of which r and d give the best result when applying Corollary 
IV.1.1.2(a, b). In case (a), n > r and r 2 > nd 2 imply n > d 2 (and even r > <i 2 ), while in 
case (b), having r 2 < nd 2 and r < n but trying to maximize r/d shows that it is enough 
to consider values of d with d < \\fn\- In short, in cases (a) and (b), we may as well 
only consider d with d < \y/n\ . Moreover, given such a d, the best choice of r is evidently 
[(iy 7 ^] for case (a) and [^V^-J f° r case ( D )- It i s still (as far as I can see) not easy to tell 
which d is best without checking each d from 1 to \y/n\ , hence I have included the script 
bestrda for case (a), and bestrdb for case (b), to do just that. Alternatively, d = [\/n\ 
often seems to be a good choice. For this choice of d and the corresponding optimal choices 
of r, (a) ends up giving a better bound than (b) if n — d 2 is even, while (b) is better if 
n — d 2 is odd. 

IV. 2. Bounds on r 

In some ways, r is easier to compute than a. For example, given Z = m(pi + ■ ■ ■ +p n ) 
for n > 9 general points, [HHF] proves by an easy specialization argument that 

r(Z)<m\V^] + \(\M-^m- 
If n > 9 is a square and m > (y/n — 2)/4, it follows (see [HHF]) in fact that 

r(Z) = m^+ r(0i-3)/2l. 

Thus r is known in some situations where a is only conjectured. 

Moreover, via an observation of Z. Ran, bounds on a give rise to bounds on r. In 
particular, given Z = m(pi + • • • + p n ) with pi general, if a(Z) > c n m for all m (where 
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c n > depends only on n), then 

T (Z) < _ 3 + \(m + l)max{v / ^, n/c}~\ 

(see Remark 5.2 of [H9]). Thus, for example, the bounds of Corollary IV.1.1.2(a, b) lead 
to bounds on r. 

It should not be surprising that r might be easier to handle than a. Being always able 
to compute ce(Z) is equivalent to being always able to compute hz and hence r(Z), while 
the reverse does not seem to be true. Moreover, arguments typically involve specializations. 
One can hope to compute r exactly using a specialization that drops a (and thereby gives 
us something to work with) while leaving r unchanged, but this of course will not work to 
compute a, only to give a lower bound. 

The scripts f indtau and unif f indtau give lower bounds for r which via the SHGH 
Conjecture are expected to be the actual values. Thus most interest is in finding upper 
bounds on r, and indeed, quite a few upper bounds have been given, both on P 2 and in 
higher dimensions (see, for example, [FL], for various results and additional references). 

Given Z = m\p\ + ■ ■ ■ + m n p n , bounding r(Z) is mostly of interest for n > 9 since 
for n < 9, for any disposition of the points, the Hilbert function of I(Z) (and hence t(Z)) 
is known (see [H4] for n < 8 or [H5]). For n > 9, the results of [H5] also allow one to 
compute r(Z) exactly, if the points pi lie on a plane cubic. If the points pi are general, 
and t is the value of r(Z') (computed via [H5]) for some specialization Z' of the points pi 
to a plane cubic, then by semicontinuity r(Z) < t. For Z = m(pi + ■ • • + p n ) with n > 9, 
this gives the bound 

t(Z) < mn/3. 

This bound is similar in concept to but better than a bound given by Segre [Seg] , obtained 
by specializing to a conic, which for Z = m{p\ + • • • + p n ) with n > 9 gives only 

t(Z) < ran/2. 

Improved bounds for r(Z) for Z = m\p\ + • • • + m n p n with pi general are given by 

Catalisano [Catl], Gimigliano [Gi3] and Hirschowitz [Hi]. For Z = ra-ipH Vm n p n with 

general points pi and mi > • • • > m n > 0, Gimigliano's result is that 

t(Z) <mH hffid 

as long as d(d + 3)/2 > n, while Hirschowitz's result is that r(Z) < d if 

\(d + 3)/2] \(d + 2)/2] > Mm + l)/2. 

i 

Catalisano's result is somewhat complicated, but generalizes and often improves 
Gimigliano's. For Z = m(pi + ■ ■ ■ + p n ) with n > 9 these all show that r(Z) is at most 
approximately m\^2n. For n sufficiently large, this clearly is better than r(Z) < mn/3. 

The bound r(Z) < m\^/ri \ + \(\^/n\ — 3)/2], mentioned above ([HHF]), results from 
specializing n > 9 points to a smooth curve of degree \y/n\ . Two other bounds which are 
also on the order of m^fn are Ballico's [B] for which r(Z) < d if 

d(d + 3) - nm(m + 1) > 2d(m - 1) - 2 
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(but note that r(Z) < m\y/n\ + \(\y/n\ — 3)/2] is better for any given n if m is large 
enough) and Xu's for which r(Z) < d if 

3(rf + 3) > (m + l)VlOn 

(although t(Z) < m\^/n\ + \(\y/n\ — 3)/2] is better if n is sufficiently large). 

By employing a sequence of specializations to infinitely near points similar to what 
he did for bounding a, Roe [R2] obtains an upper bound on r. The method applies for 
any Z = mipi + ■ • • + m n p n , with pi general and n > 2. For Z = m(pi + • ■ ■ + p n ), [R2] 
denotes this upper bound by di (m, n) and proves 

d 1 (m,n) + l < m(n/(n-l))([Tl~^((n-l + i 2 )/(n-l + i 2 -i))) < (m + l)(y/n + 1.9 + tt/8). 

The bound t(Z) < (m + l)(y/n + 1.9 + tt/8) - 1 compares very well with the bound 
r[Z) < m\\/n} + \(y/n — 3)/2] : the former is better for approximately 60% of the values 
of n between any two successive squares. 

Given a curve C, the idea of Roe's algorithm is that for any F, by taking cohomology 
of -> O x (F - C) -> O x {F) -> Cc ® 0x(F) -> 0, we have ^(I^x^)) = if 
^(X, O x {F-C)) = and ^(X, Oc^OxC-F 1 )) = 0. In Roe's case, C is always rational so 
h\X, O c ®Ox{F)) = Ois guaranteed if FC > -2, and he handles h\X, O x {F-C)) = 
by induction. 

In somewhat more detail, start with Z = m^i + • • • + m n p n with pj general, and 
so we may assume rrii > rrii+i > for all i. We have the corresponding divisor class 
F = [tEo — m\E\ — ■ ■ ■ — m n E n ] where t is as yet undetermined. Now specialize so that 
each element of S = {[Ei — E i+ i] : 1 < % < n} U {[-E n ]} and \E\ — E 2 ] is the class of a 
reduced, irreducible divisor. Now, F ■ \E\ — E 2 ] > — 1 is certainly true to start with (in 
fact, we have F ■ [E 1 — E 2 ] > 0). If F ■ [E x - E 2 - E 3 ] > -1, then fine, but otherwise 
replace F by F — [E\ — E 2 ] and unload the result with respect to S, and continue replacing 
and unloading in the same way until F ■ [E\ — E 2 — E3] > — 1. Note that throughout 
this sequence of operations we have F ■ [E\ — E 2 ] > —1, so (taking [C] = \E\ — E 2 ]) we 
have h}(X, Oc <S> Ox(F)) = 0. Also, unloading involves a succession of replacements of F 
by F — [E], where [E] is always either [Ei — E i+ \] for some i or [E n ], and can always be 
carried out in such a way that at each step we have F ■ [E] > —2. Thus we always have 
h}{E, Oe® Ox{F)) = 0, where E is a curve whose class is, at various times, [Ei — E i+1 ] 
for some % or [E n ]. 

So eventually F turns into a class for which F ■ [E\ — E 2 —E3] > —1, F • [Ei — Ei+i] > 
for all i and F ■ E n > 0. We now further specialize so that [Ei — E 2 — E3] is the class of an 
irreducible divisor, and keep replacing F by F — [Ei — E 2 — E3], unloading with respect to 
S after each replacement, as long as F ■ [Ei — E 2 — E3 — E4] < — 1. We continue in this way, 
specializing successively so that each [E\ — E 2 — ■ • • — Ej\ in turn becomes the class of an 

irreducible divisor, and replacing F by F — [Ei — • Ei] and unloading with respect to S 

after each replacement, as long as F ■ [Ei — E 2 — ■ ■ ■ — E i+ i] < — 1. Eventually we end up 

with a class F' = [tE - m' 1 E 1 m' n E n ] with F' ■ [Ei - E l+1 ] > for all i, F ■ E' n > 

and F' - [Ei- E 2 E n ] > -1. By construction, h\X, O x {F)) = for our original 

class F if h}{X, Ox(F')) = 0, but it turns out in the specialization we end up with that 
h x {X, Ox(F')) = if t > mi + m' 2 - 1. Thus Roe's bound is r(Z) < mi + m' 2 - 1. 
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By combining (in characteristic 0) the approaches of [R2] , [H9] and [HHF] , similar to 
what is done in Section IV. 1.3, Roe and I [HR] obtain another bound on r. The method 
uses a single specialization in which the same set S as above consists of classes of irreducible 
divisors, but instead of [Ei — E<i — • • • — -Ei+i] being the class of an irreducible divisor D, 
[dE — (Ei + ■ ■ ■ + E r )] is, for some d and r < n. The idea is to start with some class 
F = [tEo — m\E\ — ■ ■ ■ — m n E n ] with mi > m !+ i > for all i. We want to choose t 
to be large enough to start with so that we can keep subtracting [dEo — (E\ + ■ ■ ■ + E r )] 
and unloading with respect to S until we eventually obtain a class F' = t'[Eo] for some t' , 
while along the way always keeping /i 1 (-D, Od <S> Ox(F)) = 0. The latter is guaranteed (in 
characteristic 0) if both F ■ E > d - 2 and F ■ D > g - 1, where g = (d - l)(d - 2)/2 is 
the genus of D. 

The output of the algorithm of the previous paragraph is easy but tedious to compute 
in any given case; to get a nice formula we seem to need to choose r and d carefully. For 
example, let Z = m{pi + ■ ■ ■ + p n ) with pi general in characteristic 0. Assume r < n and 
define u > and < p < r via ran = ur + p. If 2r > n + d 2 (such as is the case for 
d = [\/n\ and r = \dy/n\), the algorithm gives 

t(Z) < max{ \(mr + g - l)/d] , (u + l)d - 2}, 

while if r < <i 2 , then the algorithm gives 

t(Z) < max{[(p + ^ - +ud, (u+l)d-2}. 

Using d = \y/n \ and r = n, the latter formula gives a bound which is always at least as good 
as that mentioned above from [HHF] . And when m is sufficiently large, the former formula 
becomes t(Z) < \mr/d+ (d — 3)/2] , which for a given n with m sufficiently large, gives a 
better bound than the bound <ii(m, n) given in [R2]. (To justify this claim, note that by a 
method similar to how [R2] shows that <ii(m, n) < — l+m(n/(n— l))(n^L 2 1 ((n — l+i 2 )/(n — 
1+i 2 — i))) one can show that m(n/(n— l))(n" = Z 2 1 ((?i— l+i 2 )/(n— 1+i 2 — i)))— Y17=3 n/(i(n— 
1) + i(i - - 2)) < d i(m,n ). But (n/(n - l^iU^dn - 1 + i 2 )/(n - 1 + i 2 - i))) > 
n/(y/n — 1 — tv/8 + l/y/n — 1); see the proof of Proposition 4.2 of [H9]. The claim now 
follows for m large enough from the fact that nj{\Jn — \ — ir/8 + l/Vn — 1) > r/d for 
d = [\/n\ and r = \dy/n\ when n > 10.) 
The formulas 

t(Z) < max{ \(mr + g - l)/d] , (u + l)d - 2} 

and 

t(Z) < max{[(p + ^ - l)/d] + ud, (u+l)d-2} 

can be evaluated by running ezunifHRtau and ezunif HRtauB, respectively. Since the 
algorithm works for any r < n and <i, it can sometimes do better than the formulas, which 
only work for certain values of r and d. Thus I have provided scripts unif HRtau and HRtau 
to compute the output of the full algorithm with respect to any specified choice of r < n 
and d. 
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V. Scripts 

We close this survey with a collection of MACAULAY 2 scripts for computing some 
of the quantities and bounds discussed above. This is a verbatim listing: There are no 
TjrjX control sequences interspersed in the text of the scripts in the TjrjXfile for this paper, 
so one can simply copy the text of the scripts from the TjrjXfile directly into a file called 
(say) BHscripts. To run a script, such as f indres (which computes a resolution of I(Z) 
for Z = m\p\ + ■ ■ • + m^Ps, where the pi are assumed to be general and each rrti is an 
integer), start MACAULAY 2 and enter the command load "BHscripts". Then enter 
the command f indres ({mi, 777.2,777.3, 777,4, 777,5,777,6,777,6,777.7,777.8}). 

The required format for each script's input parameters are described below, just before 
the listing for each script. Individual scripts can be run without loading the entire file, but 
many scripts defined below call one or more of the others, so be sure to load all scripts 
called by the one you wish to run. 

— These routines have been debugged on MACAULAY 2, version 0.8.52 

— Brian Harbourne, October 12, 2000 

— f indres: This computes the syzygy modules in any resolution 

— of the saturated homogeneous ideal defining any eight or fewer general 

— fat points of P2. The hilbert function of the ideal is also found. 

— Call it as f indres ({m_l ,m_n>) for n <= eight integers m_i . 

— Note that f indres does not rely on Grobner bases, so it is fast by comparison. 

findres = (1) -> ( 
if #1>8 then ( 

« "This script works only for up to 8 points." << endl; 
« "Please try again with an input list of at most 8 integers." « endl) 
else ( 
i:=0; 

myflag2:=0; 
w2:={>; 
ddl:=0; 
myker : =0 ; 
www : ={> ; 
ww : =1 ; 

— the list 1 of multiplicities is, for simplicity, extended if need be 

— so that it has 8 elements. 
while(#ww < 8) do ww=join(ww,{0}) ; 
n : =#ww ; 

n=n-l ; 

ww=zr(ww); — zero out negative elements of the list 

al :=f indalpha(ww) ; — find alpha, the least degree t such that I_t \ne 

dl:=al-2; 

tau:=f indtau(ww) ; 

v4:=-Q; — list of number of syzygies in each degree t listed in vO 
v3:=-Q; — list of dim of coker of \mu_t in each degree t listed in vO 
v2:={>; — list of dim of ker of \mu_t in each degree t listed in vO 
vl:={}; — list of Hilbert function values for each degree listed in vO 
v0:=O; — list of degrees from alpha-2 to tau+2, where tau is the least 

— degree such that the fat points impose independent conditions 
while (dl <= tau+2) do ( — loop from alpha-2 to tau+2, computing vO, vl and v2 

— append the current degree dl to the list of degrees 

v0=join(v0,{dl}) ; 

— append the value of the hilbert function in degree dl 
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to the list of values of the hilbert function 
vl=join(vl , {homcompdim(f undom({dl ,ww>) )}) ; 
now compute and append to the list v2 the dimension myker of the 
kernel of the map \mu_t : I_t\otimes k[P2]_l \to I_{t+1} where t=dl 
and I is the ideal of the fat points subscheme 

if dKal then (v2=join(v2 , {0}) ) ; — dKal means l_{dl}=0 so myker=0 
if dl>=al then ( — for dl>=al, compute myker 

myflag2=0; 
w2=ww ; 
ddl=dl; 

while (myflag2==0) do ( — this loop implements the main theorem of [FHH] 

— which gives an algorithm for computing myker 

w2=zr(w2) ; 
w2=prmt (w2) ; 

if homcompdim({ddl,w2})==0 then myflag2=l else ( 
if ddl*6 - ( dot (w2, {3, 2, 2, 2, 2, 2, 2,2})) <= 2 then ( 
ddl=ddl-6; 

w2 = {(w2#0)-3, (w2#l)-2, (w2#2)-2, (w2#3)-2, (w2#4)-2, 
(w2#5)-2, (w2#6)-2, (w2#7)-2}) else ( 
if ddl*5 - (dot (w2, {2, 2, 2, 2, 2, 2, 1,1})) <= 1 then ( 
ddl=ddl-5; 

w2 = {(w2#0)-2, (w2#l)-2, (w2#2)-2, (w2#3)-2, (w2#4)-2, 
(w2#5)-2, (w2#6)-l, (w2#7)-l}) else ( 
if ddl*4 - (dot (w2, {2, 2, 2, 1,1, 1,1,1})) <= 1 then ( 
ddl=ddl-4; 

w2 = {(w2#0)-2, (w2#l)-2, (w2#2)-2, (w2#3)-l, (w2#4)-l, 
(w2#5)-l, (w2#6)-l, (w2#7)-l}) else ( 
if ddl*3 - (dot (w2, {2, 1,1, 1,1, 1,1,0})) <= then ( 
ddl=ddl-3; 

w2 = {(w2#0)-2, (w2#l)-l, (w2#2)-l, (w2#3)-l, (w2#4)-l, 
(w2#5)-l, (w2#6)-l, (w2#7)}) else ( 

if ddl*2 - (dot(w2,{l,l,l,l,l,0,0,0})) <= then ( 
ddl=ddl-2; 

w2 = {(w2#0)-l, (w2#l)-l, (w2#2)-l, (w2#3)-l, (w2#4)-l, 
(w2#5) , (w2#6) , (w2#7)}) else ( 

if ddl - (dot(w2,{l,l,0,0,0,0,0,0}))< then ( 
ddl=ddl-l; 

w2 = {(w2#0)-l, (w2#l)-l, (w2#2) , (w2#3) , (w2#4) , 
(w2#5) , (w2#6) , (w2#7)}) else ( 
myflag2=2)))))))); 
if myflag2==l then myker=0 else ( 

if ddl - (dot(w2,{l,l,0,0,0,0,0,0})) == then 
myker=homcompdim({ddl-l,{(w2#0)-l, (w2#l) , (w2#2) , (w2#3) , 

(w2#4) , (w2#5) , (w2#6) , (w2#7)}})+ 
homcompdim({ddl-l,{(w2#0) , (w2#l) -1 , (w2#2) , (w2#3) , (w2#4) , (w2#5) , 
(w2#6) , (w2#7)}}) else ( 

www={3* (w2#7) +1 , 3* (w2#7) +1 , 3* (w2#7) +1 , 3* (w2#7) +1 , 
3* (w2#7) +1 , 3* (w2#7) +1 , 3* (w2#7) +1 , (w2#7) } ; 
if {8*(w2#7)+3,www}=={ddl,w2} then myker=(w2#7)+l else ( 
if homcompdim({ddl+l ,w2}) >3* (homcompdim({ddl ,w2}) ) then 
myker=0 else myker=3* (homcompdim({ddl ,w2}) ) -homcompdim({ddl+l ,w2}) ) ) ) ; 
v2= j oin ( v2 , {myker}) ) ; 
dl=dl+l) ; 

scan(#v0, i->( — this scan computes v3 from v2 and vl 
if i<2 then v3=join(v3,{0}) else ( 

if v2#(i-l)>-l then v3=join(v3 , {vl#i-3* (vl#(i-l) ) +v2#(i-l) }) else ( 
if vl#i-3*(vl#(i-l))>=0 then 

v3=join(v3,{vl#i-3*(vl#(i-l))}) else v3=join(v3 , {0}) ) ) ) ) ; 
scan(#v0, i->( — this scan computes v4 from v3 and vl 
if i<3 then v4=join(v4,{0}) else 
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v4=join(v4,{v3#i-vl#i+3*(vl#(i-l))-3*(vl#(i-2))+vl#(i-3)}))) ; 
« "The output matrix has four columns. Column 1 indicates" << endl; 
« "each degree from alpha-2 (where alpha is the least" « endl; 
« "degree t such that I_t > for the fat points ideal I)" « endl; 
« "to tau+2 (where tau is the least degree t such that the points" << endl; 
« "impose independent conditions in all degrees t or bigger)." « endl; 
« "Column 2 gives the value dim I_t of the Hilbert function" << endl; 
« "in each degree t listed in column 1. The resolution of I" « endl; 
« "is of the form -> F_l -> F_0 -> I -> 0, where F_l and" « endl; 
« "F_0 are free S=k[P2] modules. Thus F_0=oplus_t S [-t] ~{n_t}" « endl; 
« "and F_l=oplus_t S[-t]~{s_t]- for integers s_t and n_t." << endl; 
« "Columns 3 and 4 give the values of n_t and s_t in each degree t" « endl; 
« "listed in column 1 (n_t and s_t are in all other degrees)." << endl; 
transpose (matrix ({v0 ,vl ,v3 ,v4}) ) ) ) 



— findhilb: computes e(F_t(Z)) for alpha-l<= t <=tau+l, which gives 

— a lower bound for the SHGH conjectural hilbert function 

— for a fat points subscheme involving general points of P2. 

— Call it as findhilb ({m_l, ... ,m_n}) for integers m_i 

— specifying the multiplicities of the fat points in Z. 

— The conjecture is known to be correct for n<=9. 

findhilb = (1) -> ( 
ww : =1 ; 

if #1<3 then ww=join(l , {0 , , 0>) ; 

n : =#ww ; 

n=n-l ; 

ww=zr(ww) ; 

al :=f indalpha(ww) ; 

tau:=f indtau(ww) ; 

dl:=al-l; 

<< "The output gives dim I_t, computed in degrees t from alpha(I)-l to " « endl; 
<< "reg(I) , where tau = reg(I)-l is least degree such that" « endl; 
« "hilbert function of I equals hilbert polynomial of I." « endl; 
if n>9 then ( 

« "When more than 9 multiplicities are input," << endl; 
« "the output is a lower bound for dim I_t, which by the" << endl; 
« "SHGH conjecture should equal dim I_t . " « endl); 
« endl; 

« " t dim I_t" « " (tau = " « tau « ")" « endl; 

while (dl <= tau+1) do ( 

« " 11 « dl « " " « homcompdim(fundom({dl,ww>)) « endl; 

dl=dl+D) 

— input: l={n,m}, n = number of points, m = uniform multiplicity 
— output: various bounds on alpha and tau 

unifbounds = (1) -> ( 

n:=l#0; 

m:=l#l; 

ba: =bestrda(n) ; 
bb:=bestrdb(n) ; 
ea:=unif f indalpha(l) ; 
t:=0; 

<< "number of general points n of P2: " « n « endl; 
<< "multiplicity m of each point: " « m « endl; 
« endl; 

if n<= 9 then (« "Value of alpha: " « ea « endl) 
else ( 
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« "Expected value of alpha (via SHGH conjecture): " << ea << endl; 
« " Note: The SHGH conjectural value of alpha is an upper bound." << endl); 
<< "Lower Bounds on alpha:" << endl; 

<< " Roe's, via unloading: " << unif roealpha(l) << endl; 
tmp : =unif ezbhalpha(l) ; 

« " Harbourne's, via Cor IV.i.2(a, b) , using r="«tmp#l«" and d=" ; 

«tmp#2«": "« tmp#0 « endl; 

t=ea; 

while(t==unifbhalpha(l,ba#0,ba#l,t)) do t=t-l; 

<< " Harbourne's, via unloading, using r="«ba#0«" and d=" ; 

«ba#l«": "« t+1 « endl; 

tmp=ezunif HRalpha(l) ; 

« " Harbourne/Roe ' s first formula, using r="«tmp#l«" and d="« 
tmp#2«": "« tmp#0 « endl; 

<< " Harbourne/Roe ' s second formula, using r="<< (tmp#2) * (tmp#2) <<" and d="« 

tmp#2«": "«ezunifHRalphaB(l, (tmp#2) * (tmp#2) ,tmp#2) «endl; 

r:=ba#0; 

d:=ba#l; 

t=ea; 

while(t==unifHRalpha(l,r,d,t)) do t=t-l; 

t=t+l; 

tmp=ea; 

while(tmp==unifHRalpha(l,bb#0,bb#l,tmp)) do tmp=tmp-l; 

tmp=tmp+l ; 

if tmp>t then ( 

t=tmp ; 

r=bb#0; 

d=bb#l) ; 

<<" Harbourne/Roe ' s (via modified unloading), using r="«r«" and d="«d; 
« " : "«t«endl«endl; 
tt :=unif f indtau(l) ; 

if n<= 9 then ( « "Value of tau: " « tt « endl) 
else ( 

« "Expected value of tau (via SHGH conjecture): " « tt « endl; 

« " Note: The SHGH conjectural value of tau is a lower bound." « endl); 
<< "Upper Bounds on tau:" « endl; 
« " Hirschowitz's: " « Hiuniftau(l) « endl; 
« " Gimigliano ' s : " « Guniftau(l) « endl; 
if n>4 then ( « " Catalisano's: " « Cuniftau(l) « endl); 
t=0; 

while (t*(t+3)-n*m* (m+1) < 2*t*(m-l)-2) do t=t+l; 

« " Ballico's: " « t « endl; 

t=0; 

while(9*(t+3)*(t+3) <= 10*n* (m+1) * (m+1) ) do t=t+l; 
« " Xu's: " « t « endl; 
t=0; 
tmp=0 ; 

while (tmp*tmp < n) do tmp=tmp+l; 

while (2*t < 2*m*tmp + tmp - 3) do t=t+l; 

« " Harbourne/Holay/Fitchett's: " « t « endl; 

<< " Roe's, via unloading: " << unif roetau(l) « endl; 

tmp=ezunif HRtau(l) ; 

« " Harbourne/Roe ' s first formula, using r="«tmp#l«" and d="« 
tmp#2«": "« tmp#0 « endl; 

<< " Harbourne/Roe ' s second formula, using r="<< (tmp#2) * (tmp#2) <<" and d="« 
tmp#2«": "«ezunifHRtauB(l, (tmp#2) * (tmp#2) ,tmp#2) «endl; 
r=ba#0; 
d=ba#l ; 

t=unifHRtau(l,r,d,tt) ; 
tmp=unifHRtau(l,bb#0,bb#l,tt) ; 
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if tmp<t then ( 
t=tmp ; 
r=bb#0 ; 
d=bb#l) ; 

<<" Harbourne/Roe ' s (via unloading), using r="<<r<<" and d="<<d<< ": "<<t<<endl; 
if (ba#0)*(bb#0)<(ba#l)*(bb#l)*n then ( 

t = -3 + ceiling((m+l)*(bb#l)*n/(bb#0)) ; 

r=bb#0 ; 

d=bb#l) else ( 

t = -3 + ceiling((m+l)*(ba#0)/(ba#D) ; 
r=ba#0 ; 
d=ba#l) ; 

<< " Via Ran's observation, and Harbourne's bound on alpha," « endl; 
« " using r="«r«" and d="«d« ": "«t«endl) 

— input: l={ml, . . . ,mn}, n >=1 (number of points), ml, ... >=1 (the multiplicities) 

— output : various bounds on alpha and tau 

bounds = (1) -> ( 
n:=#l; 

ba: =bestrda(n) ; 

bb:=bestrdb(n) ; 

ea: =f indalpha(zr (1) ) ; 

t:=0; 

tmp : =0 ; 

<< "number of general points n of P2: " « n << endl; 
<< "multiplicities of the points: " « 1 « endl; 
« endl; 

if n<= 9 then ( « "Value of alpha: " « ea « endl) 
else ( 

« "Expected value of alpha (via SHGH conjecture): " << ea << endl; 

« " Note: The SHGH conjectural value of alpha is an upper bound." << endl); 
<< "Lower bounds on alpha:" << endl; 
« " Via Checking Psi: " « Psibound(l) « endl; 
<< " Roe's, via unloading: " << roealpha(l) << endl; 
w: =ezbhalphaA(l) ; 

« " Harbourne's, via Cor IV. i. 2(a), using r="«w#l«", d="; 
«w#2«": "«w#0«endl; 
w=ezbhalphaB(l) ; 

« " Harbourne's, via Cor IV. i. 2(b), using r="«w#l«", d="; 
«w#2«": "«w#0«endl; 
w=ezbhalphaD(l) ; 

« " Harbourne's, via Cor IV. i. 2(d), using r="«w#l«", d="«w#2; 

«", and j = "«w#3«": "«w#0«endl; 

r:=ba#0; 

d:=ba#l; 

t=ea; 

while(t==HRalpha(l,r,d,t)) do t=t-l; 

t=t+l; 

tmp=ea; 

while (tmp==HRalpha(l,bb#0,bb#l, tmp)) do tmp=tmp-l; 

tmp=tmp+l ; 

if tmp>t then ( 

t=tmp ; 

r=bb#0 ; 

d=bb#l) ; 

<<" Harbourne/Roe ' s (via modified unloading), using r="«r«" and d="«d; 
« " : "«t« endl « endl; 
tt :=f indtau(l) ; 

if n<= 9 then ( « "Value of tau: " « tt « endl) 
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else ( 

« "Expected value of tau (via SHGH conjecture): " « tt « endl; 

« " Note: The SHGH conjectural value of tau is a lower bound." « endl); 
<< "Upper bounds on tau:" « endl; 
« " Hirschowitz's: " « Hitau(l) « endl; 
« " Gimigliano ' s : " « Gtau(l) « endl; 
nn : =0 ; 
i:=0; 

w=prmt (zr (1) ) ; 

scan(#l, i->(if w#i >0 then nn=nn+l)); 

if nn>4 then ( « " Catalisano ' s : " « Ctau(l) « endl); 
<< " Roe's, via unloading: " << roetau(l) << endl; 
r=ba#0 ; 
d=ba#l ; 

t=HRtau(l,r,d,tt) ; 
tmp=HRtau(l,bb#0,bb#l,tt) ; 
if tmp<t then ( 

t=tmp ; 

r=bb#0 ; 

d=bb#l) ; 

« " Harbourne/Roe ' s (via unloading) , using r="«r«" andd="«d«": "«t«endl«endl; 
<<"0ne more bound on alpha; this one can be slow since it tries all r and d." « endl; 
w=bestbhalpha(l,ea) ; 

<< " Harbourne's alpha lower bound (via unloading), using r=" << w#l << " and d="; 
« w#2 « ": " « w#0 « endl) 

— decomp: prints a decomposition F=H+N for any divisor class F in Psi 
— as described in the print statements of the script. 
— input: decomp(l), where l={d,{ml, . . . ,mn}}, signifying the divisor 
—class F = dE_0-(m_lE_l+. . .+m_nE_n) . 

decomp = (1) -> ( 

<< "Let Psi be the subsemigroup of divisor classes generated by " « endl; 

<< "exceptional classes and by -K. For any divisor class F, this " « endl; 

<< "script determines if F is in Psi, and if so gives a decomposition" << endl; 

<< "F=H+N, where N is a sum of exceptionals orthogonal to each other and to H" « endl; 

<< "and H is in Psi but H.E >= for all exceptionals E. The point of this" << endl; 

« "is that dim |F| = dim |H|, and conjecturally dim |H| = (H~2-H.K)/2. " « endl « endl; 

i:=0; 

j:=0; 

w:=l; 

v:={}; 

ww:={}; 

ex:={}; 

mult : =0 ; 

tmp : =f undom(l) ; 

if tmp#0<tmp#l#0 or tmp#0<0 then (« "Your class is not in Psi." « endl) else ( 
if #1#1<3 then w={l#0, join(l#l , {0 , ,0» } ; 

d:=3*(#(w#l)) ; — define an element {d,v} in fundamental domain 
scan(#(w#l), i->(v=join(v,{(#(w#D) - i» ) ) ; 
ww=fundomboth(w,{d,v}) ; 

if ww#0#l == zr(ww#0#l) then (« "N = 0" « endl ) else ( 

« "N is a sum of the following fixed exceptional classes:" « endl; 

if (ww#0#0) - (ww#0#l#0) - (ww#0#l#l) < then ( 

scan(#(w#l) , j->(if j<=l then ex=j oin(ex , {1}) else ex=join(ex , {0}) ) ) ; 

ex={l , ex} ; 

ex= (f undomboth (ww#l , ex) ) #1 ; 

<<ex<<" is a fixed component of multiplicity "; 
« (ww#0#l#0) + (ww#0#l#l) - (ww#0#0) «endl ; 
ex={(ww#0#0)-(ww#0#l#l) , (ww#0#0)-(ww#0#l#0)}; 
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scan(#(w#l), j->(if j>l then ex=join(ex,{ww#0#l#j}) ) ) ; 
ww={{2*(ww#0#0)-(ww#0#l#0)-(ww#0#l#l) ,ex},ww#l» ; 
scan(#(w#l), i->( 

if (ww#0#l)#i<0 then ( 
ex={}; 

mult=-(ww#0#l)#i; 

scan(#(w#l) , j->(if j==i then ex=join(ex,{-l}) else ex=join(ex , {0}) ) ) ; 
ex={0,ex}; 

ex=(f undomboth(ww#l ,ex) )#1 ; 

<< ex << " is a fixed component of multiplicity " << mult « endl)))); 
« endl « "and H = " « (f undomboth(ww#l , {ww#0#0 ,zr (ww#0#l) }) ) #1 « endl)) 

— input: list 1 of multiplicities for fat points Z 

— output: least t such that F_t(Z) is in Psi 

Psibound = (1) -> ( 
t : =f indalpha ( 1 ) ; 
tmp:=fundom({t,l}) ; 

while (tmp#0>=tmp#l#0 and tmp#0>=0) do (t=t-l; 

tmp=f undom({t , 1}) ) ; 
t+1) 

— prat: arranges the elements of the list l={m_l , . . . ,m_n} in descending order 

— Call it as prmt(l) where 1 is a list of integers. 

prmt = (1) -> ( 
(prmtboth (1,1))#0) 

— prmtboth: arranges the elements of the list 11 in descending order, 

— and applies the same permutation to 12 

— Call it as prmt (11, 12) where 11 and 12 are lists of integers. 

prmtboth = (11,12) -> ( 

tmpvl :=11 ; 

tmpv2 : =12 ; 

vl:=ll; 

v2:=12; 

i:=0; 

j:=0; 

k:=0; 

scan(#ll, i->(scan(#ll, j->( 

if tmpvl#i < tmpvl#j then ( 
if i < j then ( 
k=-l; 
vi-O; 

v2=0; 

while (k<#ll-l) do (k=k+l; 

if k==i then (vl=join(vl , {tmpvl#j}) ; 
v2=join(v2,{tmpv2#j») else ( 
if k==j then (vl=join(vl , {tmpvl#i}) ; 

v2=join(v2,{tmpv2#i}) ) else (vl=j oin(vl , {tmpvl#k}) ; 
v2=join(v2,{tmpv2#k})))) ; 

tmpvl =vl ; 
tmpv2=v2)))))) ; 

{vl,v2}) 

— zr: replaces negative values in a list 1 by zeroes. 

— Call it as zr(l) where 1 is a list of integers. 



zr = (1) -> ( 
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v:=0; 
i:=0; 

scan(#l, i->( 

if l#i<0 then v=join(v,{0» else v=j oin(v, {l#i}) ) ) ; 

v) 

— quad: performs a quadratic transform on a divisor class dE_0- (m_lE_l+ . . . +m_nE_n) . 

— Call it as quad({d,{ml, . . . ,mn}}) . The output is 

— {2d-ml-m2-m3,{d-m2-m3,d-ml-m3,d-ml-m2,m4, . . . ,nm}}. 

quad = (1) -> ( 

i:=0; 

w:=l; 

if #1#1<3 then w={l#0 , join(l#l , {0 , , 0}) } ; 

v:={w#0 - w#l#l - w#l#2,w#0 - w#l#0 - w#l#2,w#0 - w#l#0 - w#l#l}; 
scan(#w#l, i->(if i>2 then v=j oin(v, {w#l#i}) ) ) ; 
v={2*(w#0) - w#l#0 - w#l#l - w#l#2,v>; 
v) 

— fundom: Call it as fundom({d,{ml , . . . ,mn}}) . The output is a new 

— list {d' ,{ml' mn'}}; the class dE_0-(m_lE_l+. . .+m_nE_n) is 

— equivalent via Cremona transformations to d' E_0- (m_l 'E_l+ . . . +m_n'E_n) , 

— where d' is either negative or as small as possible. 

fundom = (1) -> ( 
(fundomboth(l,l))#0) 

— fundomboth: applies fundom to 11 to reduce 11 to fundamental 

— domain of a certain group operation, and applies the same 

— group operation g to 12. If 12 starts out in the fundamental domain, 

— and {11' ,12'}=fundomboth(ll,12) , then {12 ,g~{-l}l}=f undomboth(12 ' , 1) . 

— This allows one to compute the action of g~{-l>. 

fundomboth = (11,12) -> ( 

wl:=ll; 

w2:=12; 

v:={}; 

if #11#1<3 then wl={ll#0 , join(ll#l , {0 ,0 , 0» } ; 

if #12#1<3 then w2={12#0 , join(12#l , {0 ,0 , 0» } ; 

v=prmtboth(wl#l ,w2#l) ; 

wl={wl#0,v#0}; 

w2={w2#0,v#l}; 

while ((wl#0 < wl#l#0 + wl#l#l + wl#l#2) and (wl#0 >= 0)) do ( 

wl=quad(wl) ; 

w2=quad(w2) ; 

v=prmtboth(wl#l,w2#l) ; 

wl={wl#0,v#0>; 

w2={w2#0,v#l» ; 
{wl,w2}) 

— homcompdim: computes e(F_t(Z)), the expected dimension of a component I_d 

— of a fat points ideal I corresponding to a fat point subscheme Z of general 

— points taken with multiplicities m_l, m_n. Call it as 

— homcompdim ({d, {m_l, ... ,m_n}}) ; the output is the SHGH conjectural 

— dimension of I_d, which is the actual dimension if n < 10. 

homcompdim = (1) -> ( 

h:=0; 

i:=0; 

w:=l; 
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if #1#1<3 then w={l#0, join(l#l , {0,0,0})}; 

w=fundom(w) ; 

d:=w#0; 

w=fundom({d,zr(w#l)}) ; 
d=w#0 ; 
v:=zr(w#l) ; 

if d<0 then h=0 else ( 
tmp : =0 ; 

scan(#v, i->(tmp = v#i*v#i+v#i+tmp) ) ; 
h=floor((d*d+3*d+2-tmp)/2) ; 
if h < then h=0) ; 

h) 

— findalpha: find alpha, the least degree t such that 

— I_t \ne 0, where I is the ideal corresponding to n general 

— points taken with multiplicities m_l, . . . ,m_n. Call it as 

— f indalpha({m_l , . . . ,m_n}) . The output is the SHGH conjectural value 

— of alpha, which is the actual value if n < 10 and an upper bound otherwise. 

findalpha = (1) -> ( 
i:=l; 

w:=prmt(zr(l)) ; 

if #1<3 then w=join(l, {0,0,0}) ; 

d:=w#0; — alpha is at least the max mult 

if (#w)<9 then ( — if n<=8, to speed things, make an estimate 

while (i < (#w)) do (d=d+w#i; 
i=i+l) ; 

d=ceiling(d/3) ; 

if d < w#0 then d = w#0) ; 
while (homcompdim({d,w}) < 1) do d=d+l; 
d) 

— dot: computes a dot product of two lists 11 and 12 (of equal length) 

— of integers. Call it as dot (11, 12). 

dot = (11,12) -> ( 
i:=0; 

dottot : =0 ; 

scan(#ll, i->(dottot = dottot + (ll#i) * (12#i) ) ) ; 
dottot) 

— input: l={ml, . . . ,mn}, n >=1 (number of points), ml,... (the multiplicities) 

— output: the SHGH conjectured value of tau; this is the actual value if 

— n < 10, and a lower bound otherwise. 

findtau = (1) -> ( 
t:=f indalpha(l) ; 

if t > then t = t-1; — tau is at least alpha - 1 



n 


=#1; 


V 


=1; 


P 


=dot(v,v) ; 


K 


={}; 


J 


=0; 


q 


=0; 



scan(#l, j->(q=q+l#j)) ; 

while(2*(homcompdim(join({t},{v}))) > t*t-p+3*t-q+2) do t=t+l; 
t) 

— input: positive integer n 

— output: {r,d}, where r~2>=d~2n, n>=r, and nd/r is as big as possible 



Survey of Fat Points on P 2 



bestrda = (n) -> ( 
rootn:=0; 

while (rootn*rootn<=n) do rootn=rootn+l ; 

rootn=rootn-l ; 

d:=l; 

r:=0; 

if rootn*rootn==n then r=rootn else r=rootn+l; 
tmpr : =1 ; 
tmpd: =1 ; 

while (tmpd<=rootn) do ( 
tmpr=tmpd*rootn ; 

while (tmpr*tmpr<tmpd*tmpd*n) do tmpr=tmpr+l ; 
if tmpr*d < tmpd*r then ( 

r=tmpr ; 

d=tmpd) ; 
tmpd=tmpd+l) ; 
{r,d}) 

— input: positive integer n 

— output: {r,d}, where r~2<=d~2n, n>=r, and r/d is as big as possible 

bestrdb = (n) -> ( 
rootn:=0; 

while (rootn*rootn<n) do rootn=rootn+l ; 

d:=l; 

r:=0; 

if rootn*rootn==n then r=rootn else r=rootn-l; 
tmpr : =1 ; 
tmpd: =1 ; 

while (tmpd<=rootn) do ( 
tmpr=tmpd*rootn ; 

while (tmpr*tmpr>tmpd*tmpd*n) do tmpr=tmpr-l; 
if tmpr>n then tmpr=n; 
if tmpr*d > tmpd*r then ( 

r=tmpr ; 

d=tmpd) ; 
tmpd=tmpd+l) ; 
{r,d}) 

— input: l={ml, . . . ,mn}, n= number of points, mi = multiplicity of ith point 

— output : Roe ' s algorithmic lower bound on alpha 

roealpha = (1) -> ( 

i:=0; 

v={}; 

w2:={}; 

w:=l; 

il:=2; 

if #1<3 then w=join(w, {0 , , 0» ; 
while (il<#w) do ( 
v={l>; 

scan(#w, i->(if i>0 then (if i<=il then v=join(v,{-l}) else v=join(v , {0}) ) ) ) ; 
w=zr (prmt (w) ) ; 
while (dot (w,v)<0) do ( 
w2=0; 

scan(#w, i->(w2 = join(w2 , {w#i+v#i}) ) ) ; 
w=zr (prmt (w2) ) ) ; 
il=il+l) ; 
w#0) 
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— input: l={n,m}, n >=1 (number of points), m >=1 (uniform multiplicity) 

— output : Roe ' s algorithmic lower bound on alpha 

unifroealpha = (1) -> ( 

il:=0; 

intchk: =0 ; 

n:=l#0; 

m:=l#l; 

roebnd : =m ; 

q:=n-l; — q keeps track during unloading of number of 

— points after the first with maximum multiplicity 
if n>2 then ( 
il=2; 

while (il<n) do ( 

if q<il then intchk=il*m- (il-q) else intchk=il*m; 
while (roebnd<intchk) do ( 
roebnd=roebnd+l ; 
if q<il then (q=n-il+q-l; 
m=m-l) else (q=q-il; 
if q==0 then (m=m-l; 
q=n-D); 

if q<il then intchk=il*m- (il-q) else intchk=il*m) ; 
il=il+l)) else roebnd = m; 
roebnd) 

— input: l={ml, . . . ,mn}, n= number of points, mi = multiplicity of ith point 

— output : Roe ' s algorithmic upper bound on tau 

roetau = (1) -> ( 

i:=0; 

vv={} ; 

vvl:={}; 

w2 : ={} ; 

ww : =1 ; 

il:=l; 

while (il<#ww-l) do ( 
vv={l>; 

scan(#ww, i->(if i>0 then (if i-K=il then vv=join(vv,{-l}) else vv=join(vv, {0}) ) ) ) ; 
vvl={l>; 

scan(#ww, i->(if i>0 then (if i<=il then vvl=j oin(vvl , {-1}) else vvl=join(vvl , {0}) ) ) ) ; 
ww=zr (prmt (ww) ) ; 
while ( (dot (ww,vv) ) < -1) do ( 
w2={}; 

scan(#ww, i->(w2 = join(w2,{ww#i+vvl#i}) ) ) ; 
ww=zr (prmt (w2) ) ) ; 
il=il+l) ; 
ww#0+ww#l-l) 

— input: l={n,m}-, n >=1 (number of points), m >=1 (uniform multiplicity) 

— output : Roe ' s algorithmic upper bound on tau 

unifroetau = (1) -> ( 

i:=l; 

s:=0; 

n:=l#0; 

ml:=l#l; 

m2 : =0 ; 

if n>l then m2=ml ; 

n2:=n-l; — n2 keeps track during unloading of number of 
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— points with multiplicity m2 
while (i < (n-D) do ( 

if (i+1) <= n2 then s=(i+l)*m2 else s=(i+l)*m2+n2-i-l ; 
while ((ml-s) < -1) do ( 
ml=ml+l ; 

if i < n2 then n2=n2-i else ( 
n2=n-i+n2-l; 
m2=m2-l) ; 

if (i+1) <= n2 then s=(i+l)*m2 else s=(i+l) *m2+n2-i-l) ; 
i=i+l) ; 
ml+m2-l) 

— input: 1=-Tn,m}, n >=1 (number of points), m >=1 (uniform multiplicity) 

— output: list -Ci,r,d}, with i being Harbourne's easy lower bound on 

— alpha (via Cor IV. i. 2 (a), (b)) computed using the best r and d 

unif ezbhalpha = (1) -> ( 
w:=bestrda(l#0) ; 

i:=ceiling((l#l)*(l#0)*(w#l)/(w#0)) ; — compute mnd/r rounded up 

r : =w#0 ; 

d:=w#l; 

w=bestrdb(l#0) ; 

j :=ceiling((l#l)*(w#0)/(w#D) ; — compute mr/d rounded up 
if i<j then (i=j ; 

r=w#0 ; 

d=w#l) ; 
{i,r,d» 

— input: {n,m}-, n >=1 (number of points), m >=1 (uniform multiplicity) 

— output: {a,r,d]-, where a is lower bound on alpha via formula of [HR] 

— using r and d. 

ezunifHRalpha = (1) -> ( 



n 


=1#0; 


m 


=1#1; 


t 


=0; 


d 


=0; 



while (d*d <= n) do d=d+l; 

d=d-l; 

r:=d; 

while (r*r < d*d*n) do r=r+l ; 
q:=ceiling(n*m/r) -1 ; 

while (((t+2)*(t+l)<=2*(m*n-r*q)) and t<d) do t=t+l; 
{t+q*d,r,d}) 

— input: (-fn,m},r ,d) , n >=1 (number of points), m >=1 (uniform multiplicity) 

— where r <= n, and d(d+l)/2 <= r <= d"2 

— output: lower bound on alpha via formula of [HR] . 

ezunifHRalphaB = (l,r,d) -> ( 

n:=l#0; 

m:=l#l; 

g:=(d-l)*(d-2)/2; 
tmp:=f loor((m*r+g-l)/d) ; 
t : =ceiling(n*m/r) -1 ; 
rr : =m*n-r*t ; 
s:=0; 

while (((s+l)*(s+2) <= 2*rr) and s<d) do s=s+l; 

s=s-l ; 

t=s+t*d; 
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if tmp<t then t=tmp; 
t+1) 

— input: ({n,m},r ,d) , n >=1 (number of points), m >=1 (uniform multiplicity) 

— where r <= n, and r <= d~2 

— output: upper bound on tau via formula of [HR] . 

ezunifHRtauB = (l,r,d) -> ( 

n:=l#0; 

m:=l#l; 

g:=(d-l)*(d-2)/2; 
q:=ceiling(n*m/r) -1 ; 
rr:=m*n-r*q; 

t:=q*d+ceiling((rr+g-l)/d) ; 

tmp : =q*d+d-2 ; 

if t<tmp then t=tmp; 

t) 

— input: ({n,m},r ,d,ea) , n >=1 (number of points), m >=1 (uniform multiplicity) 

— ea is an estimate for alpha (for speed) ; it must be set equal to 

— a value no bigger than the eventual value of unifHRalpha (0, for example) 

— output: lower bound on alpha via modified unloading method of [HR] , 

— computed using given r and d 

unifHRalpha = (l,r,d,ea) -> ( 
i : =ea-l ; 

n2:=0; — n2 keeps track of the number of points with maximum multiplicity 
s:=0; 
tmpi : =-1 ; 
tmpm : =0 ; 

g:=(d-l)*(d-2)/2; 
while (tmpKtmpm) do ( 

i=i+l; 

tmpi=i ; 

tmpm=l#l ; 

n2=l#0; 

if r <= n2 then s=r*tmpm else s=r*tmpm-r+n2 ; 

while (((tmpi*d-s < g) and (tmpi >= d-2)) or ( ( (tmpi+1) * (tmpi+2) <=2*s) and 

(tmpi<d) and (tmpi>=0))) do ( 

tmpi=tmpi-d; 

if r<n2 then n2=n2-r else ( 
n2=l#0-r+n2 ; 
tmpm=tmpm-l ; 
if tmpm <= then ( 
tmpm=0 ; 
n2=l#0)) ; 

if r<=n2 then s=r*tmpm else s=r*tmpm-r+n2) ) ; 

i) 

— input: ({m_l , . . . ,m_n} ,r ,d, ea) , n >=1 (number of points), m_i >=1 (multiplicities) 

— ea is an estimate for alpha (for speed) ; it must be set equal to 

— a value no bigger than the eventual value of HRalpha (0, for example) 

— output: lower bound on alpha via modified unloading method of [HR] , 

— computed using given r and d 

HRalpha = (l,r,d,ea) -> ( 

i:=ea-l; 

j:=0; 

n:=0; 

g:=(d-l)*(d-2)/2; 
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v:=prmt(zr(l)) ; 
tmpv : =v ; 
ttmpv :={}; 

scan(#l, j->(if v#j>0 then n=n+l)); 
if n>0 then ( 
ww : ={> ; 

scan(#l, j->(if j<r then ww=join(ww,{l>) else ww=j oin(ww, {0>) ) ) ; 
tmpi:=-l; 

while (tmpi < tmpv#0) do ( 
i=i+l ; 
tmpi=i ; 
tmpv=v ; 

while (((tmpi*d- (dot (ww, tmpv) )<g) and (tmpi >= d-2)) or 

(((tmpi+l)*(tmpi+2)<=2*(dot(ww,tmpv))) and (tmpKd) and (tmpi>=0))) do ( 
tmpi=tmpi-d; 
ttmpv={} ; 

scan(#l, j-> (ttmpv= j oin (ttmpv , { (tmpv# j ) - (ww# j )}))); 
tmpv=prmt (zr (ttmpv) ) ) ) ) ; 

i) 

— input: (-Cn,m}-,r ,d,ea) , n >=1 (number of points), m >=1 (uniform multiplicity) 

— ea is an estimate for alpha (for speed) ; it must be set equal to 

— a value no bigger than the eventual value of unifbhalpha (0, for example) 

— output: Harbourne's algorithmic lower bound on alpha via unloading, 

— using the given r and d. 

unifbhalpha = (l,r,d,ea) -> ( 
i:=ea-l; 

n2:=0; — n2 keeps track of the number of points with maximum multiplicity 
s:=0; 
tmpi : =-1 ; 
tmpm : =0 ; 

while (tmpi<0) do ( 
i=i+l; 
tmpi=i ; 
tmpm=l#l ; 
n2=l#0; 

if r <= n2 then s=r*tmpm else s=r*tmpm-r+n2 ; 
while ((tmpi*d-s < 0) and (tmpi >= 0)) do ( 
tmpi=tmpi-d; 

if r<n2 then n2=n2-r else ( 
n2=l#0-r+n2; 
tmpm=tmpm-l ; 
if tmpm <= then ( 

tmpm=0 ; 

n2=l#0)) ; 

if r<=n2 then s=r*tmpm else s=r*tmpm-r+n2) ) ; 

i) 

— input: ({n,m},r ,d) , n >=1 (number of points), m >=1 (uniform multiplicity) 

— where r <= n, and 2r >= n+d~2 

— output: lower bound on alpha via formula which agrees with that via unloading. 

ezunifBHalphaB = (l,r,d) -> ( 

n:=l#0; 

m:=l#l; 

q:=f loor(n*m/r) ; 

rr:=m*n-r*q; 

t :=q-l+ceiling(rr/d) ; 

tmp : =d*ceiling(n*m/r) ; 
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if tmp<t then t=tmp; 
t+1) 

— input: l={ml, . . . ,mn}, n >=1 (number of points), ml, ... >=1 (the multiplicities) 

— output: list (aa,r,d), with aa being Harbourne's lower bound on 

— alpha (via Cor IV. i. 2(a)), computed using the best r and d 

ezbhalphaA = (1) -> ( 

i:=0; 

s:=0; 

n:=0; 

v: =prmt (zr (1) ) ; 

scan(#l, i->(if v#i>0 then n=n+l)); 

w: =bestrda(n) ; 

i=0; 

while (i<n) do ( 
s=s+v#i ; 
i=i+l) ; 

best:=ceiling(s*(w#l)/(w#0)) ; 
{best, w#0, w#l}) 

— input: l={ml, . . . ,mn}, n >=1 (number of points), ml, ... >=1 (the multiplicities) 

— output: list (aa,r,d), with aa being Harbourne's lower bound on 

— alpha (via Cor IV. i. 2(b)), computed using the best r and d 

ezbhalphaB = (1) -> ( 

i:=0; 

s:=0; 

n:=0; 

v: =prmt (zr (1) ) ; 

scan(#l, i->(if v#i>0 then n=n+l)); 

w: =bestrdb(n) ; 

i=0; 

while (i<n) do ( 
s=s+v#i ; 
i=i+l) ; 

best : =ceiling(s* (w#0) / (n* (w#l) ) ) ; 
{best, w#0, w#l}) 

— input: l={ml , . . . ,mn} , n >=1 (number of points), ml, ... >=1 (the multiplicities) 

— output: list (aa,r,d,j), with aa being Harbourne's lower bound on 

— alpha (via Cor IV. i. 2(d)), computed using the best r, d and j 

ezbhalphaD = (1) -> ( 

w:={}; 

i:=0; 

scan(#l, i->(if l#i>0 then w=j oin(w, {l#i}) ) ) ; 

w=prmt (w) ; 

n : =#w ; 

bnd : =0 ; 

tmpbnd : =0 ; 

r:=0; 

d:=0; 

j:=0; 

tmpr : =0 ; 

tmpd : =0 ; 

tmp j : =0 ; 

if n>0 then ( 

while (tmpr<n) do (tmpr=tmpr+l ; 
tmpd=0 ; 
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while (tmpd*tmpd<tmpr) do (tmpd=tmpd+l ; 
tmpj=0; 

while (tmpj<tmpd*tmpd) do (tmpj=tmpj+l ; 
tmpbnd=lpa (w , tmpr , tmpd , tmp j ) ; 
if tmpbnd>bnd then ( 

bnd=tmpbnd ; 

r=tmpr ; 

d=tmpd ; 

j=tmpj))))) ; 

{bnd.r ,d, j}) 

— lpa computes the bound given in Cor IV. i. 2(d); attempts 

— various solutions with the hope of approximating the optimal 

— solution to the linear programming problem indicated by Thm IV.i.l. 

— lpa is called by ezbhalphaD 

lpa = (l,r,d,j) -> ( 

i:=0; 

n:=#l; 

sum: =0 ; 

sumb : =0 ; 

bnd:=0; 

if d*d >= r then ( 

scan(#l, i->(if i<r then sum=sum+l#i) ) ; 
bnd=ceiling(sum/d) ) else ( 
if j==0 then ( 

scan(#l, i->(if i<d*d then sum=sum+l#i) ) ; 

bnd=ceiling(sum/d) ) else ( 

M : =f loor ( (r-d*d) * (r-d*d+ j ) / j ) ; 

scan(#l, i->(if i<d*d-j then sum=sum+l#i else (if i<M+r then sumb=sumb+l#i) ) ) ; 
sumb=sumb* j / (r-d*d+ j ) ; 

if M<n-r then sumb=sumb+ (1# (M+r) ) * (r-d*d-j *M/ (r-d*d+j ) ) ; 
bnd=ceiling( (sumb+sum) /d) ) ) ; 

bnd) 

— input: ({ml, . . . ,mn},r,d,ea) , n >=1 (number of points), ml, ... >=1 

— (the multiplicities), r and d positive integers, ea any value 

— not bigger than the eventual value of bhalpha; can be set to 

— output: Harbourne's unloading lower bound on alpha, using given r and d 

bhalpha = (l,r,d,ea) -> ( 

i:=ea-l; 

j:=0; 

n:=0; 

v : =prmt (zr (1) ) ; 
tmpv : =v ; 
ttmpv:={}; 

scan(#l, j->(if v#j>0 then n=n+l)); 
if n>0 then ( 
ww:={>; 

scan(#l, j->(if j<r then ww=join(ww,{l>) else ww=j oin(ww, {0>) ) ) ; 
tmpi:=-l; 

while (tmpi < tmpv#0) do ( 
i=i+l ; 
tmpi=i ; 
tmpv=v ; 

while ((tmpi*d- (dot (ww, tmpv) )<0) and (tmpi >= tmpv#0)) do ( 
tmpi=tmpi-d; 
ttmpv=0 ; 

scan(#l, j-> (ttmpv=join(ttmpv, { (tmpv#j ) - (ww#j ) }) ) ) ; 
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tmpv=prmt (zr (ttmpv) ) ) ) ) ; 

i) 

— Find bhalpha using best possible r and d; 

— ea is an a priori estimate for alpha (for speed) 

— it must be set to a value >= than the actual value 

— of alpha (e.g., ea=f indalpha(l) ) 

bestbhalpha = (l,ea) -> ( 

w:={}; 

i:=0; 

scan(#l, i->(if l#i>0 then w=join(w, {l#i}) ) ) ; 

w=prmt (w) ; 

n : =#w ; 

bnd : =0 ; 

tmpbnd : =0 ; 

r:=0; 

d:=0; 

tmpr : =0 ; 

tmpd : =0 ; 

if n>0 then ( 

while (tmpr<n) do (tmpr=tmpr+l ; 
tmpd=0 ; 

while (tmpd*tmpd<tmpr) do (tmpd=tmpd+l ; 
tmpbnd=ea; 

while (tmpbnd==bhalpha(w, tmpr, tmpd, tmpbnd) ) do tmpbnd=tmpbnd-l ; 

tmpbnd=tmpbnd+l ; 

if tmpbnd>bnd then ( 

bnd=tmpbnd ; 

r=tmpr ; 

d=tmpd) ) ) ) ; 

{bnd,r,d}) 

— input: ({ml , . . . ,mn} ,r ,d, et) , n >=1 (number of points), ml, ... >=1 (the multiplicities), 

— r and d positive integers, et any lower bound for tau (used for speed; can be set to 0) . 

— output: Harbourne/Roe ' s algorithmic upper bound on tau, 

— with given r and d (assumes char = 0) . 

HRtau = (l,r,d,et) -> ( 

i:=et-l; 

j:=0; 

n:=0; 

v:=prmt(zr(l)) ; 
tmpv : =v ; 
ttmpv : ={} ; 

scan(#l, j->(if v#j>0 then n=n+l)); 
if n>0 then ( 
ww:=-Q; 

g:=(d-l)*(d-2)/2; — genus of plane curve of degree d 

scan(#l, j->(if j<r then ww=join(ww,{l]-) else ww=j oin(ww, {0}) ) ) ; 

tmpi : =0 ; 

while (tmpv#0 > 0) do ( 
i=i+l; 
tmpi=i ; 
tmpv=v ; 

while ((tmpi*d- (dot (ww, tmpv) )>=g-l) and (tmpi>=d-2) and (tmpv#0 >0)) do ( 
tmpi=tmpi-d; 
ttmpv=-Q ; 

scan (#1 , j -> (ttmpv= j oin (ttmpv , { (tmpv# j ) - (ww# j )}))); 
tmpv=prmt (zr (ttmpv) ) ) ) ) ; 
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i) 

— input: (-fn,m},r ,d,et) , n >=1 (number of points), m >=1 (the uniform multiplicity), 

— r and d positive integers, et any lower bound for tau (used for speed; can be set to 0) 

— output: Harbourne/Roe ' s algorithmic upper bound on tau, using given r and d 

— (assumes char = 0) . 

unifHRtau = (l,r,d,et) -> ( 
i:=et-l; 

n2:=0; — n2 is the number of points with maximum multiplicity 
s:=0; 
tmpm: =1 ; 
tmpi : =0 ; 

g:=(d-l)*(d-2)/2; — genus of plane curve of degree d 
while (tmpm > 0) do ( 

i=i+l; 

tmpi=i ; 

tmpm=l#l ; 

n2= 1#0; 

if r<=n2 then s=r*tmpm else s=r*tmpm-r+n2 ; 
while ( (tmpi*d-s>=g-l) and (tmpi>=d-2) and (tmpm >0)) do ( 
tmpi=tmpi-d; 

if r<n2 then n2=n2-r else ( 
n2=l#0-r+n2; 
tmpm=tmpm-l ; 
if tmpm<0 then ( 

tmpm=0 ; 

n2= 1#0)) ; 

if r<=n2 then s=r*tmpm else s=r*tmpm-r+n2) ) ; 

i) 

— input: l={n,m}, n >=1 (number of points), m >=1 (the uniform multiplicity) 

— output: list {a,r,d>, where a is Harbourne/Roe ' s formulaic upper 

— bound on tau (char 0) computed using r and d 

ezunifHRtau = (1) -> ( 

n:=l#0; 

m:=l#l; 

d:=0; 

while (d*d <= n) do d=d+l; 

d=d-l; 

r:=d; 

while (r*r < d*d*n) do r=r+l ; 

g:=(d-2)*(d-l)/2; 

a:= ceiling( (m*r+g-l)/d) ; 

b:=-2+d*ceiling(m*n/r) ; 

if a<b then a=b; 

{a,r,d}) 

— input: l={n,m}-, n >=1 (number of points), m >=1 (the uniform multiplicity) 

— output: the SHGH conjectured value of alpha; this is the actual value if 

— n < 10, and an upper bound otherwise. 

unif f indalpha = (1) -> ( 

n:=l#0; 

m:=l#l; 

a:=-l; 

if n==l then a=m; 
if n==2 then a=m; 
if n==3 then a=ceiling(3*m/2) ; 
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if n==4 then a=2*m; 

if n==5 then a=2*m; 

if n==6 then a=ceiling(12*m/5) ; 

if n==7 then a=ceiling(21*m/8) ; 

if n==8 then a=ceiling(48*m/17) ; 

if n==9 then a=3*m; 

if n>9 then ( 

while (a*a-n*m*m+3*a-n*m+2 <0) do a=a+m; 

a=a-m ; 

while (a*a-n*m*m+3*a-n*m+2 <=0) do a=a+l) ; 

a) 

— input: l={n,m}, n >=1 (number of points), m >=1 (the uniform multiplicity) 

— output: the SHGH conjectured value of tau; this is the actual value if 

— n < 10, and a lower bound otherwise. 

uniffindtau = (1) -> ( 



n: 


=1#0; 






m: 


=1#1; 






t: 


-i; 






if 


n==l 


then 


t=m-l ; 


if 


n==2 


then 


t=2*m-l; 


if 


n==3 


then 


t=2*m-l; 


if 


n==4 


then 


t=2*m; 


if 


n==5 


then 


t=ceiling((5*m-l)/2) ; 


if 


n==6 


then 


t=ceiling((5*m-l)/2) ; 


if 


n==7 


then 


t=ceiling((8*m-l)/3) ; 


if 


n==8 


then 


t=ceiling((17*m-l)/6) ; 


if 


n==9 


then 


t=3*m; 


if 


n>9 then 


( 



while (t*t-n*m*m+3*t-n*m+2 <0) do t=t+m; 
t=t-m; 

while (t*t-n*m*m+3*t-n*m+2 <0) do t=t+l) ; 
if t<0 then t=0; 
t) 

— input: l={n,m}, n >=1 (number of points), m (the multiplicity of each point) 

— output: Hirschowitz ' s lower bound for tau 

Hiuniftau = (1) -> ( 

n:=l#0; 

m:=l#l; 

t:=m; 

s:=n*m* (m+1) ; 
a:=ceiling((t+3)/2) ; 
b:=ceiling((t+2)/2) ; 
while (a*b*2 <= s) do ( 
t=t+l; 

a=ceiling((t+3)/2) ; 
b=ceiling((t+2)/2)) ; 

t) 

— input: l={ml, . . . ,mn}, n >=1 (number of points), ml, ... >=1 (the multiplicities) 

— output: Hirschowitz ' s lower bound for tau 

Hitau = (1) -> ( 

n:=#l; 

i:=0; 

w: =prmt (zr (1) ) ; 
t : =w#0 ; 
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s:=0; 

scan(#l, i->(s=s+(w#i)*((w#i)+l))) ; 
a:=ceiling((t+3)/2) ; 
b:=ceiling((t+2)/2) ; 
while (a*b*2 <= s) do ( 
t=t+l; 

a=ceiling((t+3)/2) ; 
b=ceiling((t+2)/2)) ; 

t) 

— input: l={n,m}, n >=1 (number of points), m (the multiplicity of each point) 

— output: Gimigliano's lower bound for tau 

Guniftau = (1) -> ( 

n:=l#0; 

m:=l#l; 

t:=0; 

while (t*(t+3)<2*n) do t=t+l; 
m*t) 

— input: l={ml, . . . ,mn}, n >=1 (number of points), ml, ... >=1 (the multiplicities) 

— output: Gimigliano's lower bound for tau 

Gtau = (1) -> ( 
n:=0; 

w:=prmt(zr(l)) ; 

scan(#l, i->(if w#i >0 then n=n+l)); 

t:=0; 

s:=0; 

i:=0; 

while (t*(t+3)<2*n) do t=t+l; 
scan(#l, i->(if i<t then s=s+w#i)); 
s) 

— input: l={n,m}-, n >=5 (number of points), m>0 (the multiplicity of each point) 

— output: Catalisano's lower bound for tau 

Cuniftau = (1) -> ( 

s:=l#0; 

m:=l#l; 

r:=0; 

t:=0; 

f :=0; 

while(f*(f+l) <= 2*s) do f=f+l; 
f=f-l; 

while (2*r<2*s-f*(f+l)) do r=r+l; 

dl:=0; 

d:=f ; 

if r==0 then dl=f-l else dl=f ; 
t=dl+(m-l)*d; 

if 2*t+l < 5*m then t=ceiling( (5*m-l) /2) ; 
if t<2*m-l then t=2*m-l; 

if r == f then (if s >= 9 then t=m*dl+l) ; 
t) 

— input: l={ml, . . . ,mn}, n >=5 (number of points), ml, ... >=1 (the multiplicities) 

— output: Catalisano's lower bound for tau 

Ctau = (1) -> ( 
n:=0; 
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i:=0; 

w: =prmt (zr (1) ) ; 

scan(#l, i->(if w#i >0 then n=n+l)); 

vm:={}; 

vs:=-Q; 

i=0; 

while (i < n-1) do ( 

if w#i > w#(i+l) then ( 

vs=join(vs , {i+1}) ; 

vm= j oin ( vm , {w#i}) ) ; 

i=i+l); 
vs=join(vs,{n>) ; 
vm=join(vm, {w#(n-l)}) ; 
i=#vm - 1 ; 
v: =-[vm#i]- ; 
while (i > 0) do ( 

v=join({vm#(i-l) - vm#i},v); 

i=i-l); 
vf :={}; 
vr:={}; 

scan(#vs, i->( 
f :=0; 
r:=0; 

while (f*(f+l) <= 2*(vs#i)) do f=f+l; 
f=f-l; 

while(2*r<2*(vs#i)-f*(f+l)) do r=r+l; 
vf =join(vf ,{f }) ; 
vr= j oin ( vr , {r» ) ) ; 
t:=0; 

if (vr#(#vr-l)) == then t = - 1; 
dl:=t+vf#(#vf-l) ; 
scan(#vs, i->( 

t=t+(vf#i)*(v#i))) ; 
if 2*t+l < (w#0)+(w#l)+(w#2)+(w#3)+(w#4) then 
t=ceiling( ( (w#0) + (w#l) + (w#2) + (w#3) + (w#4) -1) /2) ; 
if t<(w#0)+(w#l)-l then t=(w#0) +(w#l) -1 ; 

if (vr#0) == (vf#0) then (if s >= 9 then (if (w#0)==(w#(n-l) ) 

then (if (w#0)>l then t= (w#0) *dl+l) ) ) ; 
if (vr#0) == then (if s > 9 then (if (w#0)==(w#(n-2) ) 

then (if (w#(n-l)) == 1 then t=(w#0) *dl+l) ) ) ; 

t) 
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